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The lepton flavor violating interactions are worthwhile to examine since they 
are sensitive to physics beyond the Standard Model. The simplest extension of 
the Standard Model promoting the lepton flavor violating interactions are the 
so called two Higgs doublet model which contains an additional Higgs doublet 
carrying the same quantum numbers as the first one. In this model, the lepton 
flavor violating interactions are induced by new scalar Higgs bosons, scalar 
and pseudo scalar A^, and Yukawa couplings, appearing as free parameters, are 
determined by using the experimental data. On the other hand, the possible 
extra dimensions are interesting in the sense that they ensure a solution to the 
hierarchy and cosmological constant problems and also result in the enhancement 
in the physical quantities of various processes. In the present work, we predict 
the branching ratios of lepton flavor violating radion decays r e^,fi^, r 
e^,r^ and r —> /i^,r^ in the two Higgs doublet model, including a single extra 
dimension, in the framework of the Randall Sundrum scenario. We observed that 
the branching ratios of the processes we study are at most at the order of 10"*^ 
for the small values of radion mass and it decreases with the increasing values of 
the radion mass. Among the LFV decays we study, the r — /i^, decay would 
be the most suitable one to measure its branching ratio. 
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CHAPTER 1 



INTRODUCTION 



The major goal of physics has always been the simplification and the unification of 
seemingly diverse and complicated natural phenomena. In the second half of the 
twentieth century, as a result of successful approaches, a significant progress has 
been made in the particle physics, in the identification of fundamental particles 
and the unification of their interactions. Glashow-Weinberg-Salam [1, 2] com- 
bined the quantum electrodynamics (QED) and the weak interactions into the 
electroweak (EW) theory and the Standard Model (SM) of elementary particles 
has emerged, which can be considered as a good example satisfying this major 
goal of physics. Being a quantum field theory (QFT) based on the gauge group 
SU{3)c ® SU{2)l ® U{1)y, the SM describes all of the known elementary con- 
stituents of the universe together with the three out of four fundamental forces: 
the strong force, the weak force, and the electromagnetic force. In the QFT, all 
interactions are mediated by means of force carrier particles, mediators. In the 
case of electromagnetic interaction the mediator is the photon (7), one of the 
four gauge bosons of the group SU{2)l ^ f/(l)y, and for the weak interactions, 
there exist the remaining gauge bosons of the same group, the so called W^, 
bosons. In addition to this, for the strong interactions, the mediators are eight 
gluons (Gj), the gauge bosons of the group SU{3)c- The remaining force, called 
as the gravity, is far too weak to be of any consequence at the experimentally 
accessible energy scales that are relevant to particle physics. In addition to the 
mediator particles, the SM contains matter particles: the Higgs boson and the 
fermions, namely leptons and quarks which fall into three generations. The first 
generation contains all stable stuff of which the stable matter is composed. The 
second and third generations of particles decay, therefore, they are not present 
in the stable matter and physicists are still trying to understand their role in the 
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underlying theory. For each quark and lepton there exists a corresponding anti- 
quark and antilepton. This is all adding up to an embarrassingly large number 
of elementary particles: 12 leptons, 36 quarks, 12 mediators, and, as we will see 
later, Glashow-Weinberg-Salam theory calls for at least one Higgs particle, so we 
have a minimum of 61 particles to content with. In the next chapter, we will 
see how this structure leads to the first consistent and self-contained theory. The 
energy range which defines this theory extends up to several hundreds GeV. For 
the details of the model construction see for example textbooks [3, 4], and the 
review [5] existing in the literature. 

The SM has been very successful in explaining many diverse experimental 
results in the energy range available at present. However, behind this energy 
range it possesses some conceptual problems which motivate us to look physics 
beyond. The big issues in physics beyond the SM can be conveniently grouped 
into four categories. The Unification: What is the reason beyond the hierarchy 
of fundamental forces? The problem of Flavor: Why are there so many different 
types of quarks and leptons? The Mass problem: What is the origin of masses of 
fundamental particles and their mass hierarchies? Does the Higgs boson exist? 
The cosmological constant problem. In addition to the conceptual problems of 
SM, there also exist phenomenological hints obtained from measurements of fla- 
vor changing neutral currents (FCNCs), including lepton flavor violating (LFV) 
interactions which also indicate the need for physics beyond the SM since the SM 
predictions differ from the upper limits coming from current measurements. 

There are various alternative extended models proposed for solving these prob- 
lems of the SM such as the multi Higgs doublet model (MHDM) [6, 7, 8, 9], 
the minimal supersymmetric model (MSSM) [10, 11, 12, 13], left-right (super) 
symmetric model [14], the Zee Model [15], the see-saw model [16], technicolor 
model [17], extra dimensional models [18]-[31]; large extra dimensions [19, 20, 21], 
universal extra dimensions (UED) [22, 23, 24], non-universal extra dimensions 
(NUED) [25, 26], split fermion scenarios [27, 28, 29], the Randall- Sundrum model 
(RS model) [30, 31]. 

Based on the phenomenological hints, the violation of flavor symmetry in 
the leptonic sector is of special interest to physicists. In the SM, the FCNCs 
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with massless neutrino, are not allowed in the lepton sector and, in the quark 
sector, they are prohibited at tree level, despite they seem not to violate any 
fundamental law of nature. The negligibly small branching ratios (BRs) of the 
decays based on the FCNCs stimulate one to go beyond the SM and they are 
worthwhile to examine since they open a window to test new models, to ensure 
considerable information about the restrictions of the free parameters, with the 
help of the possible accurate measurements. An elegant framework to open up 
the possibility of the tree level FCNCs is proposed through the general two Higgs 
doublet model (2HDM) (see [7, 8, 9] for details), the most primitive candidate of 
MHDM, which is obtained by adding a second Higgs doublet, having the same 
quantum numbers as the first one. This doublet may lead to FCNCs in its Yukawa 
sector, representing interactions between the Higgs fields and fermions (see for 
example [9]). In this model, the lepton flavor (LF) violation is driven by the 
new scalar. In addition, the mass hierarchy problem among third generation of 
quarks, namely the top and bottom quarks, also could be solved in the scope of 
2HDM, such that, unlike in the SM where both quarks gain mass through the 
interaction with the same Higgs doublet, there is a possibility that the bottom 
receives its mass from one doublet (say 0i) and the top from the another one (say 
02)- Then the hierarchy of their Yukawa couplings could be more natural. 

A theory which consists of the SM, combined with gravity, contains two enor- 
mously different energy scales. One is the EW scale niEw ~ 10^ GeV at which 
EW symmetry is broken, and the other is the Planck scale Mpi ~ 10^^ GeV 
which determines the strength of gravitational interactions. Newton's laws state 
that the strength is inversely proportional to the second power of that energy, 
and because the strength of gravity is so small, the Planck scale mass (related 
to the Planck scale energy hj E = mc^) should be very large. Generally, when 
making predictions in particle physics, we can ignore gravity since the gravita- 
tional effects on particles in the EW energy scale are completely negligible. But 
that is precisely a question which particle physicists try to find an answer: Why 
is the gravity so weak? A solution to this problem comes from models with ex- 
tra dimensions where gravity becomes strong and cannot be neglected. In 1998, 
Nima Arkani-Hamed, Savas Dimopoulos and Gia Dvali [20, 21] proposed a model 
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(called as ADD Model) with n compact extra spatial dimensions of large size to 
bring the Planck scale down to TeV scale. Depending on the details of their im- 
plementation, the space in their model contains two, three or more compact extra 
dimensions. For two extra dimensions, the hierarchy problem in the fundamental 
scales could be solved and the true scale of quantum gravity would be no more 
the Planck scale but of the order of EW scale. This is the case that the gravity 
is spreading over all the volume including the extra dimensions and thus it is 
diluted by a large volume of them so much that it will be very feeble in the lower 
dimensional effective theory^, although it is very strong in higher dimensions. On 
the other hand, the matter fields together with the electromagnetic, strong and 
weak forces are restricted in four dimensions, called four dimensional (4D) brane. 
Unlike gravitational force, these forces will not be accessible to the higher dimen- 
sions. As mentioned above, in ADD model, the extra dimensions are compact 
and their compactification leads to the appearance of towers of heavy Kaluza- 
Klein (KK) modes [32] of particles such that, in 4D effective theory, the existence 
of the extra dimensions are felt by the appearance of these KK modes. How- 
ever, since the matter fields do not travel along the extra dimensions but bound 
to 4D brane, they will not carry extra dimensional momenta. In other words, 
none of the SM particles will have the KK partners. The only particle that will 
have KK partners is the graviton, the force carrier particle of the gravitational 
force. Since the KK partners of graviton also interact with gravitational strength 
(i.e., as weakly as graviton itself), it would be no easier to produce or detect KK 
partners of the graviton than to observe the graviton itself which also has never 
directly seen by anyone up to now. This means that, gravity, being the only force 
which lives in higher dimensions, the existence of large extra dimensions will not 
contradict with the experimental results. Despite the success of ADD proposal 
in solving the hierarchy problem, there exist also some weaknesses of the theory. 
In fact their model do not actually solve the hierarchy problem, because one still 
have to solve why the size of extra dimensions are so large. 

An alternative approach is introduced by Randall and Sundrum (RSI model) 

^Effective theory is a theory describing those elements and forces that are in principle ob- 
servable at the distance or energy scales over which it is applied. 
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[30, 31] to explain the huge discrepancy between ruEw and Mpi without the need 
for a large extra dimension, or for any arbitrary large number at all. In this 
scenario, the geometry is a non-factorizable one where the gravity is localized 
in a 4D brane, so called Planck brane, which is one of the boundary of the 
extra dimension and away from another 4D brane, TeV brane, which is the other 
boundary where we live^. Theory also includes a finely tuned 5D cosmological 
constant A which serve as sources for 5D gravity and in 4D, with the help of 
opposite tensions on boundaries, it vanishes. 

The review topics we include in this thesis is, broadly, divided into three cat- 
egories: In Chapter 2, we give a brief review of the SM. Chapter 3 is devoted to 
the simplest extension of the SM, the so called the 2HDM. In Chapter 4, we give 
a summary of models with extra dimensions. In Chapter 5, we investigate the 
branching ratios of lepton flavor violating radion decays r — > e^,/i^, r e^,r^ 
and r — >■ h^,t^ in the 2HDM, in the framework of the Randall Sundrum sce- 
nario (RSI). Chapter 6 represents our conclusions. In Appendix A, we present 
the global and local gauge invariance. Appendix B is devoted to the detailed cal- 
culations of Einstein equations, that we use in our review. Appendix C represents 
the calculation of spin connection. 



^The extra dimension is compactificd to /Z2 orbifold with two 4D brane boundaries which 
reside at the orbifold fixed points. 
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CHAPTER 2 



THE STANDARD MODEL 

In the second half of the twentieth century physicists made an impressive contri- 
bution to the progresses in particle physics with complementary theoretical and 
experimental studies, whereupon the SM of elementary particles has emerged 
[1, 2]. The SM, being a QFT (see for example [33]), describes all of the known 
elementary particles together with the three out of four fundamental interac- 
tions of nature. According to the SM, the elementary particles constituting the 
universe are called as fermions (i.e., they have spin one-half), namely quarks 
and leptons and the fundamental interactions are the electromagnetic force, the 
weak force (responsible for radioactive decay) and the strong force (which holds 
atomic nuclei together). The SM is based on the principle of gauge symmetry 
(see Apendix A), which means that the properties and interactions of elementary 
particles are governed by certain symmetries which are related to the conserva- 
tion laws. Therefore, the electromagnetic, weak, and strong forces are all gauge 
forces and they are mediated by the exchange of certain particles, called gauge 
bosons (i.e., they have spin one) which are the photon, the and bosons, 
and eight gluons, respectively. Several attempts have been made to fit gravity, 
the remaining force, into this gauge framework but these attempts are resulted in 
failure. However, the gravity is too weak to change particle physics predictions 
in the current experimental energy scales. 

With these in mind, it is worthwhile summarizing the Fermi theory [34] which 
describes the weak interaction phenomenology in the mid-1950. 
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2.1 The Fermi Theory 



The progress in the field theory of the weak interaction was rather stagnant for 
many decades, from Fermi's attempt to describe the (3 decay 

n ^ p + + Ue, (2.1) 

in 1933, to the advent of the gauge theories in the 1970s. Fermi expressed this 
decay mathematically as, at a single point in the space-time, the quantum me- 
chanical wavefunction of a neutron is transformed into the wavefunction of a 
proton, and that the wavefunction of an incoming neutrino is transformed into 
that of an electron. He wrote the phenomenological Lagrangian as 

CF = ^jlU^)jr,^,{x) + h.c., (2.2) 

where Gp = (1.16639 ± 0.00002) x 10"^ GeV~^ is the Fermi coupling constant, 
J^Md^^) = iJp{x)T ai)n{x) and J^^pti^) = ^e(a;)r"^i,^(x). This action, from very 
start, was known to suffer from a series of problems. First of all, the Fq, ma- 
trices, that contain the essence of the weak interaction, consist of all possible 
combinations of the 16 Dirac matrices. It took many years to narrow down the 
choice. In 1958, Feynman and Gell-Mann [35] with the help of further experi- 
mental data proposed that the correct combination of Fq, matrices should only 
contain a mixture of vector and axial- vector^ (V-A) quantities written in the form 
Tq = 7a (1^75) to incorporate the parity- violating effects of the weak interaction. 

Since the weak force is of extremely short range, Fermi's theory of point-like 
interaction yields excellent approximate results at low energies. However, at high 
energies the theory suffers from additional problems. The main problem is the 
violation of unitarity. Uf, + e~ ^ Uf, + e~~ scattering is one of the simplest example 
of the weak interaction processes. The Feynman diagram for this scattering in 
the Fermi's picture of four point interaction is given in the Fig. 2.1. In the 

^an axial vector (or a pseudovector ) is a quantity that transforms like a vector under a 
proper rotation, but gains an additional sign flip under an improper rotation (a transformation 
that can be expressed as an inversion followed by a proper rotation). 
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Figure 2.1: The four point + e — ^ i^e + e scattering. 

center of mass frame (CM), the differential cross section is found as 

da _ Glk"^ 

with CM four momentum k, and 

AG^k^ 



(2.3) 



(2.4) 



TT 

where k^ ^ m^. Since the four fermion interaction takes place at a single point in 
space-time, the differential cross section is a pure s-wave. Partial wave unitarity 
for s-wave requires that 

Using the above equation, k is obtained as 

Then, above a certain energy (i.e., k > 300 GeV), Fermi theory violates unitarity. 
Thus, we can say that it is an effective theory up to the energies k < 300 GeV. 

Another problem in this theory is the non-renormalizability. Unfortunately, 
even in the lowest order approximation in four fermion interactions, one encoun- 
ters horrible divergences which cannot be eliminated by proper renormalization. 
The Feynman diagram of the four point interaction including the lowest order 
correction for the scattering z/g + e~ ^ z/g + e~ is shown in the Fig. 2.2 below 
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which is generated by multiplying the current-current interaction with itself and 
the amplitude for this scattering is proportional to oc 4^ = oo"^. To eliminate 




Figure 2.2: One loop correction to the Ue + e — > z/g + e scattering . 

these problems, the idea is that the weak interaction is mediated by intermediate 
massive vector boson exchange. The Feynman diagram for this process is shown 
in the Fig. 2.3. Therefore, we replace Cp defined in the eq. 2.2 with 




Figure 2.3: The + e ^ Uf. + e scattering mediated by intermediate massive 
vector boson exchange. 

= gwr{x)Waix) + h.c, (2.7) 

where Wa{x) is the weak intermediate vector boson. Now, in the lowest order 
diagram the differential cross section of the u^ + e" — > z/g + e" scattering is found 
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as 

da_ ^ 2g^e 

for again k"^ > m^. Here, m^y and q are the mass and the momentum transfer 
vector of the W boson, respectively. As 0, the new Lagrangian in the eq. 

2.7 reduces to Fermi Lagrangian given in the eq. 2.2 provided 

However, the interaction is no longer point-like, but mediated by the force carrier 
particles. In the scope of the SM, all fundamental interactions are mediated by 
the exchange of gauge bosons. These interactions together with the corresponding 
gauge bosons which mediate the forces are listed in the following table in order 
of decreasing strength. 



Table 2.1: The four fundamental forces in nature. 



Force 


Strength 


Range 


Theory 


Mediator 


Strong 


10 


< lO^i^m 


Chromodynamics 


Gluon 


Electromagnetic 


10-2 


oo 


Electro dynamics 


Photon 


Weak 


10-13 


< lO^i^m 


Flavordynamics 


W^,Z 


Gravitational 


10-42 


oo 


Geomet r o dynamics 


Graviton 



According to the QFT, the short range of the weak force could mean only one 
thing: the weak gauge bosons had to have non zero masses. The mechanism that 
gives rise to the masses of gauge bosons is known as the Higgs mechanism [36] 
which relies on the phenomenon of spontaneous symmetry breaking (SSB) which 
we will consider in the following section. 



2.2 Spontaneously Broken Symmetries 

SjTumetry is one of the most important aspects of theoretical particle physics, 
since the basis of our current description of nature originate in symmetries so 
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that every continuous symmetry leads to a conservation law. A system is said 
to be symmetric if it remains invariant after applying a set of transformation 
rules that constitute a mathematical group. Symmetries are categorized into two: 
spatial symmetry and internal symmetry. In the case of spatial symmetry physics 
threats all directions and all positions as the same, internal symmetries tell us 
that physical laws act the same way on distinct, but effectively indistinguishable 
objects. The fundamental forces, electromagnetic, weak, and strong forces all 
involve internal symmetries. (Gravity is related to the symmetries of space and 
time). Exact symmetries are fairly rare in nature. Thus, the symmetries that the 
usual 4D theories possess can be broken explicitly or spontaneously. In the case 
of SSB of gauge symmetries, which is one of the crucial ingredients of the SM, if 
the broken symmetry is global, the Goldstone theorem [37] applies, whereas if it 
is local, then we have Higgs mechanism [36]. In general, the phenomenon of SSB 
is simply stated as follows. 

"A system is said to possess a symmetry that is spontaneously broken if the 
ground state of a dynamical system does not possess the same symmetry properties 
as the Lagrangian". Here, the ground state -its vacuum state- is the state in which 
the field has its lowest possible energy. Now, we will use a toy model (see [4] for 
details) to explain the Goldstone theorem and the Higgs mechanism. 

2.2.1 The Goldstone theorem 

Let us consider the following Lagrangian density which describes a couple of self 
interacting complex scalar fields (j){x) = (pi{x) + i4>2{x) and its complex conjugate 
(f)*{x) = (f)l{x) - i(f)*2{x) 

£ = {d,<p){dy*) - - A(00*)^ (2.10) 

where //^ is regarded as the bare mass of the field quanta and A is the term 
for self interaction. It is clear that the Lagrangian remains invariant under the 
group U{1) of global gauge transformations (see Appendix A for details). For 
SSB, we should check whether the ground state of the system will be invariant 
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under global gauge transformations or not. For constant (j) the kinetic term, 
{d^(t)){d^(f)*) vanishes. Then, the ground state is obtained when the potential term 
y (0, (jf) = fi'^(j)(f)* + A(00*)^ corresponds to the minimum. Since the potential 
term is a function of (p and </>* only in the combination of (fxp*, we can make a 
change of variables so that, p = (pcj)*. Substituting this into V{(f), 0*) we get 

Vip)=p'p + Xp\ (2.11) 



The minimum of the potential can be obtained only if A > 0, which we take to 
be so. However, p"^ can have both positive and negative values if we do not insist 
on interpreting p as mass. To find the minimum of the potential, we take the 
derivative of V{p) with respect to p and equate this derivative to zero such that. 



dVjp) 
dp 



p^ + 2Ap = 0. 



(2.12) 



Since p = (fxj)* = (0i + z02)(0* — «02) — il'Pil'^ + I'/'sP), P can only take positive 
values. Therefore, for p'^ > a. unique minimum occurs at the origin p = 0, 
i.e., at = and the function V{p) looks like as in figure 2.4. Then, we have 

V(rho) 



rho 



Figure 2.4: The potential function for positive /i^ 



a symmetric ground state configuration under the group U{1) of global gauge 
transformations for /i^ > 0. On the other hand, for p^ < 0, = is not a 
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minimum. Instead, the minimum is at p = — ;U^/2A, i.e., at |</)| = v / \/2 with 
V = \/ —[i^ I A. Any value of satisfying this relation will give us a true ground 
state such that 

= ^e^^, (2.13) 

where A is real. Then, we have a continuum degenerate set of ground states for 
negative values of yU^. In this case, the function 1^(0) looks like as in figure 2.5. 
Each will not be symmetric under the global gauge transformation defined in the 



V(phi) 




Re(phi) 



Figure 2.5: The potential function for negative /i^. 

eq. A-2. Then, using the definition of SSB made above one can simply conclude 
that the symmetry of the Lagrangian has been spontaneously broken. We are 
free to choose any point on the ring of minima since they are equivalent. Let us 
choose this point to be on the real axis such that 

= + (2.14) 

where ^(x), x{^) ^-^e real fields and ^(x) = x(a;) = in the ground state. Substi- 
tuting into the eq. 2.10 and ignoring constant terms we get 

= \mf + \{d'xf - \v\' - \vi{e + x') - \Ke + x^f. (2.15) 
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Then, we end up with a massless x{^) fi^ld and a field ^{x) with a spontaneously 
generated mass 

m^{x) = 2Xv'^. (2.16) 



Let us now examine a less trivial example [5] given by 



(2.17) 



where (p is an n-component real scalar field and C is invariant under the orthogonal 
group in n dimensions, 0{n). Again for /i^ < we find a whole ring of minima 
whenever 0*0* = —fi'^/X is satisfied. In this case, we are free to choose one of 
the 0* to be non-zero in the ground state. Let it be the n'^^ component of such 
that, 

/o\ 





(2.18) 





\v J 

The number of generators that original symmetry group 0(n) possesses is ^n{n — 
1). There is also a non-trivial subgroup 0(n — 1), which has ^{n — l){n—2) number 



of generators leave the vacuum invariant. Let Lj,- be the ^ 



n(n 



1) independent 



matrices that generates 0{n) and hjlkj = Lij for i,j ^ n] be the ^{n — l){n — 2) 
matrices generating 0{n — 1). There remains — 1 independent matrices which 
are denoted by ki[ki = L^n] with 1 < i < n — 1. Defining r] and with again 
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l<z<n — Iwe get 



/ \ 





Note that in general, 



For j =n we have 



\ V + r] J 



{Lij)kl — —i[^ik^jl — ^il^jk\- 



{Lin)kl — {ki)kl — —ilSikSnl — ^ii^nfc]- 

Operating ki on the column vector Vi = v6in we get 



{ki)jivi 

-i[SijSnl - 5il6nj]Vi 
-i[6ijVn - SnjVi] 
-i[6ijV6nn - SnjVSin] 



-iv6., 



(2.19) 



(2.20) 



:2.2r 



(2.22) 



Thus, in the lowest order 0j = < n) and (j)n = v + rj. Then, the Lagrangian 
density in terms of and t) can be presented as 

^ = \[d^r^d,v) + d^i,d,i,] - ^fi'iv + - X{v + v)' + ... (2.23) 

Looking at this Lagrangian, we can say that the field t] has a positive mass of 
— 2/i^ and the (n — 1) scalar fields remain massless. These massless bosons are 
called as Goldstone bosons. In conclusion, for every broken generator that leaves 
the vacuum invariant there exists a massless Goldstone boson. 
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2.3 The Higgs Mechanism 



Now, we will use the same Lagrangian in the eq. 2.10 but impose invariance 
under U{1) of local gauge transformations (see Appendix A for details). To make 
the Lagrangian invariant under this transformation we must replace the partial 
derivative by the covariant derivative = — ieA^ and add a kinetic term 
— jF^yF^^. Then, the Lagrangian density becomes 

£ = -\f^,,F^^'' + [{d^ + teA^^)ct>*{d, - teA^)<P] - - K<P*4>)\ (2-24) 

where is the massless gauge field. Under local gauge transformations we have, 

0(x) ^0'(x) =exp-^^(^)0(x), 

0*(a;)_,0*'(a;)=exp^^(-)0(x), (2.25) 
A,-.A^{x)=A,{x)-U,e{x). 

Again we will look at the minimum in the potential. For A > and > we 
obtain a symmetric ground state at </> = 0. However, when fi^ < there exists 
again a ring of degenerate ground states. Proceeding as before we set 

^{x) = ^[v + ^{x)+zx{x)], (2.26) 

with V = \/ —^^ I A so that (\)vac = v/V^- Substituting this into 2.24, we obtain 
the following Lagrangian density 

4 z z z (2 27) 

-l{2Xv^)i^-evA,d^X + - 

It is surprising that the gauge field A^ seems to acquire mass in the quantum 
picture. The Lagrangian density in the eq. 2.27 now seems to describe the 
interaction of a massive gauge field A^ and two scalar fields. To ensure the gauge 
invariance, the gauge transformations in terms of ^(x) and x{x) should be in the 
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following form 



^ (x) ^'{x) = [v + ^ (x)] COS 9{x) + x{x) sin 9{x) — v, 

Xix) x'ix) = x{x) cos e{x) -[v + ^(x)] sin 6'(x), (2.28) 

A,^ A'^{x)=A,{x)--^d,e{x). 

We are free to choose 6{x) to be the phase of (f){x) since the theory is invariant 
under any choice of transformation of this function. Then, 

= exp-'^(") = ^[v + r]{x)], (2.29) 
V2 

will be real, with ri[x) is real. Substituting these into the eq. 2.24 the Lagrangian 
density becomes 

C = - -Fl^F'^'^ + \d^7^d,7^ + \e'v^A'^A>^ 

(2.30) 

with F'^j^F'^^ = dfj^A'j^ — d^A'^- By writing the Lagrangian density in this form, 
we can say that it describes the interaction of the massive vector field A'^ with 
the massive, real, scalar field rj. This field is called as Higgs field with a mass 
of 2Xv^ = -2/^2. In this way, all massless particles completely disappears. Con- 
sequently, in spontaneously broken symmetries the gauge boson acquires mass 
due to disappearing Goldstone boson. Therefore, for each massive gauge fields 
we need a complex scalar field, one piece of which disappears and reappears as 
the longitudinal mode of the vector field. Scalar part of this complex field, the 
so called Higgs boson, remains. 

2.4 The Standard Model Lagrangian 

Having discussed the ingredients of the SM, let us turn our attention to the SM 
Lagrangian, Csm- As mentioned above the SM is based on the principle of gauge 
symmetry. The overall gauge group of the SM, under which the SM Lagrangian 
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remains invariant, contains both the Quantum Chromo dynamics (QCD) and the 
unified EW interaction and is written symbohcally as 



GsM = SU{3)c ® SU{2)l ® U{1)y. 



(2.31) 



The first group, SU{3)c, represents QCD. The subscript C indicates that the 
gauge bosons of QCD, the eight gluons, couple only to colored particles, quarks. 
The remaining part SU{2)l x f/(l)y represents the EW interaction, proposed 
by Glashow-Weinberg-Salam [1, 2], with the subscripts L and Y indicating that 
the group SU{2)l couples only left handed particles and that the group f/(l)y 
couples to weak hypercharged particles where the hypercharge is obtained using 
the Gell-Mann-Nishijima relation [38] Q = T3 + Y/2. The EW theory, developed 
by Glashow-Weinberg-Salam, predicted the masses of the gauge bosons W"^ and 
Z° to be about 80 GeV and 90 GeV, respectively. In 1983, physicists at CERN 
[39] led by Carlo Rubia were able to produce and measure the masses of the 
and which were in complete agreement with the predictions of EW theory. The 
discovery of these particles may be considered as the first experimental evidence 
of the SSB. In the minimum formulation of the SM, a complex scalar doublet 
(the Higgs field) is required, which is denoted by 0, so that by interacting with 
the gauge bosons, it produces the desired breaking 



where U{1)q is a subgroup of S'f/(2)i®f/(l)y. This breaking of symmetry occurs 
due to the non-zero vacuum expectation values (VEV) of the scalar field of the 
form 



The reason for why this type of breaking occurs is as follows. The and Z° 
gauge bosons are massive. Therefore, SU{2)l ^ U{1)y can not be a symmetry 
of the vacuum, whereas the photon, being massless, refiects that U{1)q is a good 
symmetry of the vacuum. 



SU{3)c^U{1)q, 



(2.32) 




(2.33) 
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Now, let us write the most general renormalizable EW SM Lagrangian. It can 
be divided into five parts: 



C-SM ^kinetic ~^ ^kinetic ~^ ^kinetic ~^ ^pot ~^ ^ ■ (2.34) 

The CJ^-in^iif. term corresponds to the fermionic sector of the SM Lagrangian. 
It includes both the left-handed and right-handed chiralities and can be presented 
as 

>CLet.c = E ^Ltl^'D.i^L + ^R^^'D.i^n, (2.35) 

where ipL = |(1 ~ Ts)"^ stands for the left-handed weak isodoublets and t/^r = 
|(1 + 75)^/^ for the right-handed weak isosinglets. The normal derivative, 9^, is 
replaced by the covariant derivative, D^: 

D, = d^ + igWln + i^-!^B^Y, (2.36) 

to preserve the gauge invariance. Here, g and g' are the coupling constants 
associated with the groups SU{2)l and U{1)y gauge groups, respectively. The 
corresponding generators to each gauge group are and Y in order. Moreover, 

are the three weak interaction bosons and is the single hypercharge boson. 
Here, the gauge boson fields W^, W^, couple to weak isospin and couple 
to weak hypercharge. 

The second part of the SM Lagrangian is the kinetic term for the scalar Higgs 
field, and is responsible for the interaction of the gauge and Higgs fields. It can 
be written as 

^LeHc=pM0)^P/.0), (2.37) 

with is defined in the eq. 2.36. 

The corresponding kinetic term for the gauge fields reads 

lA„„..„. 1 



^kinetic = ~4 PrPliv - -^B^'^B^,^, (2.38) 



i=l 
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where = df_iW^ — d^W^^ — ge^^^WjyV^ is the antisymmetric field strength 
tensor of the group SU{2)l with e*-''^ being the group structure constant and 
Bfj,u = dp^By — dpB^ is that of the group t/(l)y. After a proper normahzation of 
the gauge fields, the photon, the neutral weak boson, and the charged weak 
boson W"^ fields are obtained as 

= sinOwW^ + cosOwBf,, 

= cos 9wWj^ - sin 9wB^, (2.39) 



where 



sin 6w = — , ^ ; cosOw = , ^ =, (2.40) 



with 9w being the weak mixing angle. Finally, the photon becomes massless and 
the mass eigenstates for and Z^ bosons are obtained as 



M^.=f ; M,.. = ^t±ll. (2.41) 
The Higgs potential denoted by Cf^'^ in the SM Lagrangian reads 

<o* = /i'(0V) + A(0V)', (2.42) 

where /i^ and A are the free parameters. For /i^ < the scalar field develops a 
non-zero vacuum expectation value at |0|„ac = v/\/2 with v = \/ —y? jX. Thus, 
the Lagrangian gains a set of ground states for negative values of ^ . As a result, 
symmetry of the Lagrangian is spontaneously broken. In addition, through the 



Higgs mechanism, the Higgs mass is yielded to be equal to mn = y/2Xv. Notice 
that, all the terms explained up to now in the SM Lagrangian are CP invariant. 

The final piece of the SM Lagrangian, £y*^ called the Yukawa Lagrangian 
describes the interaction among the fermions and the Higgs field. The general 
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form can be expressed as^ 



Cy = vFA^L<pDj^R + vl'jQi,L4>Uj,R + vfjkL<pEj,R + h.c, (2.43) 

where 4> = ^^20* and r]fj^'^^s are responsible for the masses of up-down quarks 
and leptons, respectively. In addition, Qi^i, Uj^R and Dj ^ denote the left handed 
doublet, right handed up and right handed down quarks, respectively. Similarly, 
h,L represent the left handed leptons and ^ the right handed ones. These 
fermions are presented in a more elegant way in the following table: 



Table 2.2: The known fermions. 



Generation 


Quarks 


Leptons 




Charge 2/3 


Charge —1/3 


Charge —1 


Charge 




Color (R G B) 


Color (R G B) 


Colorless 


Colorless 


First 

Mass(GeV) 


u (up) 

0.0015 -0.003 


d (down) 
0.003 - 0.007 


e (electron) 
0.000511 


z/g (electron neutrino) 
< 3 X 10-9 


Second 
Mass(GeV) 


c (charmed) 
1.25 ±0.09 


s (strange) 
0.095 ±0.025 


/i (muon) 
0.106 


(muon neutrino) 
< 190 X 10-6 


Third 

Mass(GeV) 


t (top) 
174.2 ±3.3 


b (beauty) 
4.2 ±0.07 


r (tau) 
1.777 


Vt- (tau neutrino) 
< 18.2 X 10-3 



In this table, quark masses given correspond to the approximate rest mass 
energy of quarks confined in hadrons since free quarks have not been observed yet. 
For each quark and lepton given in the table there is a corresponding antiquark 
and antilepton. 

As mentioned before, in the SM, the fundamental fermionic constituents of the 
matter are quarks and leptons. All properties of these particles are summarized 
in the Table 2.2. These particles are placed into SM as left-handed doublets and 

^In the case of massive neutrinos, there is an additional term in the Lagrangian: = 
riijliXfiiVj^R + h.c. 
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right-handed singlets. The left handed doublets for quarks are given by 




(2.44) 



and the right handed singlets for quarks are 

dR ; Ur ; Sr ; Cr ] bn ] tR. (2.45) 
On the other hand, the lepton doublets are 

(:). ^ (:;). ■ (:•).■ 

and the right handed singlets for leptons are 

eR ] /^R ; Tr. (2.47) 

In the charged weak interactions of leptons, the coupling of takes place 
strictly within a particular generation in the case of massless neutrinos . In 
other words, upper members of left handed lepton doublets couple to the lower 
members in the same doublet. That is, only the vertices e~VeW~ , ^~v^W~, 
and t~Ut-W~ appear , however, there is no cross generational vertices such as 
e~v^W~ . The coupling of to quarks is not quite so simple since there exist 
cross generational vertices as well, such as suW~ . The idea is that, the quark 
generations are rotated for the purposes of weak interactions such that 




(2.48) 



where d! , s', and 6', the linear combinations of the d, s, and b, are obtained by 
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using the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix Vij 
( d'\ I v.... K,. \ ( d\ 



s' 
b' 



Vud Vus Vub 

Vcd Vcs Vcb 

Vtd Vts Vtb 



(2.49) 



where the off-diagonal elements of the CKM matrix allow flavor transitions be- 
tween different generations. The experimentally measured values of the matrix 
elements are [40] 



/ 0.97377 ±0.00027 0.2257 ± 0.0021 0.00431 ± 0.0003 \ 
0.230 ±0.011 0.957 ±0.017 ±0.093 0.0416 ± 0.0006 
0.0074 ±0.0008 0.0406 ±0.0027 > 0.78 



(2.50) 
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CHAPTER 3 

BEYOND THE STANDARD 

MODEL 

The SM has been extremely successful in describing the behavior of all known 
particles in elementary particle physics up to the EW energy scale at the order 
of 10^ GeV. However, behind this energy scale, it possesses some conceptual 
problems which motivates us to look physics beyond the SM. There are various 
alternative extended models proposed for solving these problems as indicated 
in the introduction part. MHDM [6] is one of them. In this chapter, we will 
introduce the simplest extension of the SM, the so called the 2HDM [7, 8, 9]. 

3.1 The Two Higgs Doublet Model 

Let us first present the motivation for examining the 2HDM; 

• In the SM, it is assumed that the Higgs sector must be minimal having only 
one physical neutral Higgs scalar. However, there is no fundamental reason 
favoring this minimal choice. The 2HDM, being the simple extension of the 
SM, possesses five physical Higgs bosons, namely, a charged pair (i^^), two 
neutral CP even scalars {H^ and /i"), and a neutral CP odd scalar (A^). 

• The ratio between the masses of top and bottom quarks is rrit/mi, ^ 174/5 ~ 
35. According to the SM, both quarks gain mass through interactions with 
the same Higgs doublet. Then, we end up with an unnatural hierarchy 
between the corresponding Yukawa couplings. In the scope of the 2HDM, 
there is a possibility that the bottom receives its mass from one doublet 
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(say (pi) and the top from the another one (say (j)2)- Then the hierarchy of 
their Yukawa couphngs would be more naturaL 



• In the framework of the SM the flavor is conserved in the lepton sector 
for massless neutrinos. The LFV interactions, carried by the FCNCs, exist 
in the extended SM, the so called uSM, at least at one loop level, which 
is constructed by taking the neutrinos massive and permitting the lepton 
mixing mechanism [41, 42]. However, even in the z/SM, due to the smallness 
of the neutrino masses, the theoretical predictions of the BRs of the LFV 
interactions are too small to reach the experimental limits. In addition, in 
the quark sector the FCNCs are prohibited at tree level despite they seem 
not to violate any fundamental law of nature. In that aspect, the 2HDM 
is an elegant framework to open up the possibility of the tree level FCNCs 
in both lepton and quark sectors which is driven by the new scalar Higgs 
bosons S, the CP even scalar h^, and the CP odd scalar A^, and controlled 
by the Yukawa couplings. 

• The 2HDM is a minimal extension in that it adds the fewest new arbitrary 
parameters. Instead of one free parameter of the SM, this model has six 
free parameters: the four Higgs masses, the ratio of the VEVs, tan/5, and 
a Higgs mixing angle, a. Notice that vf + f| is fixed by the W mass 

Having stated our motivations for an additional scalar doublet, our next task is 
to introduce the 2HDM. In the 2HDM, a second Higgs doublet having the same 
quantum numbers as the first one is introduced such that. 
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with hj^ercharges Y = 1. Parameterizing the doublets in a more convenient way 
we can write them in the following form 

^1 = I ) ; ^2 = ( I , (3.2) 

V2 / \ V2 



with the VEVs 




('^•i) = ^ ; ('^2) = ^ , (3.3) 




where v = {vf + v^y/^ = {V2Gf)~^^^ = 246 GeV. Here, and are the 
CP even, is the CP odd neutral Higgs bosons, and is the charged Higgs 
boson. 

In the 2HDM, the Higgs part of the SM Lagrangian should be extended to 
include the interaction with the second Higgs doublet. Then, the kinetic term in 
eq. 2.37 becomes 

(D^$i)t(D^$i) + (D^<|.2)t(D'^$2), (3.4) 

and the most general renormalizable CP invariant Higgs potential potential is 
written in the form 

y($i, $2) = Ai($I$i - vlf + A2($^$2 - vjf 

+ A3[($I$i-t;?) + (<^>5<^'2-t^2')F 

+ A4[($I$i)($^$2) - ($l$2)($^$i)] (3.5) 

+A6[/m($^$2)]^ 

where the parameters Aj are real. 

Then, what remains is the Yukawa piece of the SM Lagrangian in the presence 
of the two scalar doublets which is written as follows: 

(3.6) 

+ ri^h^L^iEj^R + ^ffi^L^2Ej,R + h.c, 
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where for i = 1,2 are the two scalar Higgs doublets, $j = ^o"2$^, rjfj^'^ 
and i^-^'^ are off diagonal 3x3 matrices of the Yukawa couplings where i, j 
denote family indices (see Chapter 2 for the definitions of terms appearing in the 
Lagrangian). 

As stated before, FCNCs with massless neutrino are naturally suppressed in 
the tree level in the SM. To avoid FCNCs at tree level, one can explicitly impose 
the following ad hoc discrete symmetry sets 

(/) $1 ^ $2 ^ $2, Dj^R -Dj,R, Uj^R -Uj^R, 

{II) $1 ^ -$1, $2 ^ $2, Dj^R ^ -Dj- iJ, Uj,R +Uj,R, 

(3.7) 

into the Yukawa Lagrangian, Imposing these discrete symmetry sets, 

the so called model I and model II are obtained depending on whether the up- 
type and the down-type quarks are coupled to the same or two different scalar 
doublets, respectively. If no discrete symmetry is implemented into the C^hdmi 
both up-type and down-type quarks and also leptons will have flavor changing 
(PC) couplings. This type of 2HDM is called as model III where we should take 
into account all the terms in the Yukawa Lagrangian given in eq. 3.6. 

It is possible to make a rotation of the doublets in such a way that only one 
of the doublets acquire VEV so that 

j ; < $2 >= 0. (3.8) 
The two doublets in this case arise of the form 

V2 / \ V2 
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Here H'^ and are not the neutral mass eigenstates. The neutral mass eigen- 
states are obtained from {H^, H^, H"^) as follows 



H = [{H^ — v)cosa — sina], 

= [-{H^ -v)sina + cos a], 
A° = H'^, (3.10) 

where a is the mixing angle. It is also possible to express and as 

functions of mass eigenstates 

= {Mucosa — h^sina) + v, 
= {h^cosa + sina), 

= v4°. (3.11) 

Choosing a = 0, ifi becomes the well known mass eigenstate h^. As mentioned 
before, the model III version of the 2HDM opens up the possibility of FCNCs at 
the tree level and the FC part of the Yukawa Lagrangian reads 

£"0"" = egg,L$2f/,,ii + ^SQ,,L^2D,,n + if~kL^2E,,R + h.c. (3.12) 

There is another version of the 2HDM, the so called Model IV such that 0i 
couples and give masses to up-type quarks and 02 couples and gives masses to 
the down-type quarks and the VEV of the doublets are chosen 

< $1 >= I ° j ; < $2 >= (| J^^^ 1 . (3.13) 

In this case the term containing Ag is replaced by XQ{I'm{(f)f(f)2) —V1V2 sin^)^) and 
is responsible for the CP violation in the Higgs sector. 



28 



CHAPTER 4 



EXTRA DIMENSIONS 



In recent years, models with extra dimensions, not yet experienced and not yet 
entirely understood, have been studied extensively in the literature (see for ex- 
ample [18]- [31]). The strong motivation to study these scenarios comes from 
the fact that they resolve some of the most basic mysteries of our universe such 
as the hierarchy problem between the two fundamental energy scales, the EW 
scale {rriEw ~ 10'^ GeV) and the Planck scale {Mpi ~ 10^^) GeV, where the 
strength of gravity becomes comparable to the one of other gauge interactions. 
In fact, there is an important difference between these two energy scales. While 
the electroweak interactions have been probed at distances ~ m~^^ gravitational 
interactions has not remotely been probed at distances ~ Mpl: it has only been 
accurately measured in the 1 cm range. Apart from the hierarchy problem, the 
cosmological constant problem (for reviews see [43, 44]), the puzzle of why the 
vacuum energy is driven to be a very small number, has also been tried to be 
solved within the extra dimensional scenarios. One possible explanation for the 
smallness of the cosmological constant problem can be stated as, if there is a 4D 
theory with only 4D sources, these will necessarily lead to an expanding universe. 
However, if there is 4D sources in 5D, the effects of brane sources can be balanced 
by a 5D cosmological constant to get a theory where the effective 4D cosmolog- 
ical constant would be vanishing. In this way, for an observer on a brane, the 
universe will still appear to be static and flat. The 5D background, however, will 
be curved since there exists a bulk cosmological constant in there. Such extra 
dimensional theories are called as warped extra dimensions, where the brane is 
kept fiat and the extra dimensions are curved. This solution to the cosmological 
constant problem first pointed out by Rubakov and Shaposhnikov [45]. Finally, 
the extra dimensional scenario, named as the split fermion scenario [27, 28, 29], 
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provides an alternative view for the fermion mass hierarchies by assuming that 
the fermions were located at different points in the extra dimensions with the 
exponentially small overlaps of their wavef unctions. 

Among the models with extra dimensions, emergence of the ordinary four 
dimensional SM as the low energy effective theory of more fundamental theory 
lying in higher dimensions found acceptance. The process of passing from a 
fundamental theory to the effective theory includes the compactification of the 
extra dimension(s). This compactification leads to the appearance of towers of 
heavy KK [32] modes of particles in 4D effective theory. 

4.1 Large Extra Dimensions 

In 1998, Nima Arkani-Hamed, Savas Dimopoulos and Gia Dvali [20, 21] tried to 
explain the hierarchy problem between the scales niEw and Mpi by assuming 
the existence of n extra compact spatial dimensions of large size. According to 
this model (called as ADD Model), the SM fields are confined to the 4D brane 
while gravity is free to propagate in large extra dimensions. In other words, the 
gravitational field has extra components in n large extra dimensions and this 
extra components cause it to be weaker than the other forces at long distances 
because it would have been diluted by the large volume of the extra dimension. 
In this scenario, since the SM particles are confined to a 4D brane, everything 
that does not involve gravity would look exactly the same as it would without 
the extra dimensions, even if the extra dimensions were extremely large. These 
extra dimensions, being compact, lead to the appearance of towers of heavy KK 
modes. However, since the SM particles, which are confined to a brane, would 
not have KK partners. The only particle in this model that must have the KK 
partner is the graviton. However, the graviton's KK partners interact far more 
weakly than the SM KK partners (see for example [47] and the references therein). 
Therefore, the KK partners of graviton would be much harder to be observed. 
One question ADD wanted to address with their set up is how large the size of the 
extra dimensions can be without getting into conflict with observations made up 
to date. To answer this question we need to match the 4D effective theory to the 
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fundamental higher dimensional one. Let us assume that the higher dimensional 
action takes the same form with Einstein-Hilbert action: 



Si+n ~ j rf^+"xv/^(4^i?(^+"). (4.1) 



Here \/g^^ and are the metric^ and curvature scalars in 4 + dimen- 

sions. The action should be dimensionless and we need to multiply the higher 
dimensional Lagrangian by the appropriate power of the fundamental Planck 
scale (M^,) , to make the action dimensionless. Since terms \J g(^~^^) and 

carry dimensions of —n — 4, and 2, respectively, the power of has to 
be n + 2. Thus, we write the higher dimensional action as 

= -M;+' J d^+"xv/^(^i?(^+"). (4.2) 

Now, we will try to find how the usual 4D action 

^4 = -Ml, J (fx^/^)R^'^\ (4.3) 

is contained in the higher dimensional expression. For time being, it is assumed 
that the space-time is flat, and the extra dimensions are compact. So the metric 
is written as 

ds'^ = ivt^u + h^u)dx''dx'' - r'^dnf^^, (4.4) 

where four dimensional coordinate with /i = 0,1,2,3, dQ"^^^ is the line 

element of the fiat extra dimensional space, r corresponds to the radius of the 
extra dimension, ?7^j, is the fiat 4D metric and h^^i, is the 4D fluctuation of the 
metric around its minimum. Since our goal is to find out how a 4D theory emerges 
from a higher dimensional one, we put only 4D fluctuation into the metric. Then, 



^Metric is a quantity that establishes the measurement scale that determines the physical 
distances and the angles. A metric on 3D space can take the form ds^ = axdx^ +aydy^ +azdz'^, 
where x, y. z are the three coordinates of space, and ax, ay, az can be numbers or functions 
of X, y, z. If ax = ay = = 1, we have flat space. More complicated metrics can have cross 
terms, such as dxdy. In that case, the metric must be described by a tensor with two indices 
that is denoted by the coefficients which is the coefficient in front of dxidxj. 
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the expressions for ^/g^^ and are obtained as 

Substituting these expressions into the eq. 4.2 we get 

= -M:+2 y y d'x^/^R^'\ (4.6) 

where the factor J dVti^n)'^'^ is simply the volume of the extra dimensional space, 
V(n). Comparing the eqs. 4.3 and 4.6, we obtain the matching relation for the 
gravitational couplings as 

Mil = M:+V(„) ~ r"M:+2. (4.7) 

Now, we will find the matching relation for gauge couplings. Based on the as- 
sumption that every field propagates in all dimensions we can write the action 
as 

= - ! d^'^'^x^FMNF''''^/^^). (4.8) 
J 4:gi 

Here, Fmn with M,N = 0,l,...3-|-?7,isthe higher dimensional field strength 
tensor and is the fundamental gauge coupling. Taking the integral over the 
extra dimension we get 



^(4) 



I d'x^F,,F^^^). (4.9) 



Then, we obtain the matching of gauge couplings as 

1 K 



glff 9t 



(4.10) 



Notice that, in the eq. 4.8, d^~^"'x, ^/g^^^ and Fmn carry dimensions —n — 4, 
and 2, respectively. Thus, the mass dimension of g^, has to be —n/2 so that the 
action remains dimensionless. Looking at the mass dimensions of g^ and we 
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can write 



9* n. (4.11^ 



Substituting this into the eq. 4.10 we get 

1 



2 -KM:~r"M:. (4.12) 

Seff 



Using the eq. 4.7 we obtain an equation for M^, as 



2 

11^ n + 2 

M, ~ (4.13) 

n+2 



If we substitute this into the eq. 4.12 we get 



2n 

1 M^+' 



f,2 



(4.14) 



Then, r becomes 

1 2i+2 

r - ■ (4.15) 

Since r ~ 1/Mpi there would be no hope of finding out about the existence of 
these tiny extra dimensions in the near future. 

Up to now, we have assumed that every field propagate in all dimensions. 
Now, we will deal with the restrictions on the size of extra dimensions when the 
SM particles are localized to 4D brane while the unobserved fields such as gravity 
are to propagate in extra dimensions. However, in distances as small as the size of 
extra dimensions, it is impossible to test gravity. Because, in such short distances 
electromagnetic and weak forces become more dominate than the gravitational 
force and the gravitational interactions have been tested at the distances of the 
order of millimeter. Therefore, the real bound on the size of extra dimensions 
becomes 

r < 0.1mm, (4.16) 

if only gravity propagates in the extra dimension. From the relation Mp; ~ 
]\^n+2^n know that if we increase the value of the radius r, the fundamental 
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Planck scale M,, decreases. If < ITeV, we would have observed some effect of 
quantum gravity in the collider experiments. Thus, one has to impose the lowest 
possible value to be M* ~ ITeV. Therefore, one can say that being equal to the 
fundamental Planck scale, niEw is the only fundamental short distance scale in 
nature and Mpi is valid in a 4D scenario is an effective scale derived from niEw- 
Such models are called as Large extra dimensions, proposed by Arkani- Hamed, 
Dimopoulos and Dvali. 

Let us check, how large a radius one would need for the lowest possible value 
of M*, using the eq. 4.7 the value of radius would be 

- = M,(^)^ = (ITel-)lO-". (4.17) 
r Mpi 

Using 

lGeV~^ = (4.18) 

we get 

r ~ 2.10"^^10~^cm. (4.19) 

For n = 1, the value of r = 2.10^^ cm. Since this value is larger then the 
astronomical unit of 1.5 x 10^^ cm, we can conclude that there cannot be one 
flat large extra dimension. If there are two extra dimensions, r ~ 2 mm. This is 
just a borderline for the currently probed gravitational experiments. For n > 2, 
r < 10~^ cm. This value is so small to be measured in the near future. Then, the 
hierarchy problem between the fundamental Planck scale and the scale of weak 
interactions would have been solved so that gravity would be weaker than the 
other forces at long distances because it would have been diluted by the large 
volume of the extra dimension. 
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4.2 Universal and Non-Universal Extra Dimen- 
sions 



As mentioned above, the compactification of the extra dimensions to a circle 5*1 
with a small radius R makes them imperceptibly small. But can an extra dimen- 
sional universe hide its nature so completely that none of its features distinguishes 
it from a 4D world? That would be hard to believe. If there are extra dimen- 
sions, fingerprints of them sure to exists. Such fingerprints are called as KK [32] 
particles. These new particles originate in extra dimensions, but appear to us 
as extra particles in our 4D space-time. In other words, they are manifestations 
of particles, which are in higher dimensions, in 4D such that every particle that 
travel in higher dimensional space is replaced by KK particles in our 4D space. 
A universe with extra dimensions contains both familiar particles and their KK 
relatives that carry extra dimensional momentum. However, a 4D dimensional 
space-time does not include information about higher dimensional position or 
momentum. This extra dimensional momentum would be seen in our 4D world 
as mass. Thus, KK particles should be like the ones we know (having the same 
charge), but heavier. If the universe contains additional dimensions, these heavier 
KK particles will be the first real evidence of them. 

The wavefunction of a KK particle is written as Fourier decomposition of the 
higher dimensional wavefunction. To be concrete let us imagine a space with 
only one additional spatial dimension which is compactified on S^/Z2 orbifold 
{Z2 : y —y)- Then, the KK decompositions of the five dimensional (5D) fields 
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can be written as follows, 



n=l 



AJx,y) = + cos(^) 



n=l 



(4.20) 



n=l 



f/(^,y) = 4^4°^ (^) + ^f:[f^!r^(^)cos(^) + t/f)(x)sin(|^)], 

* * n=l 

= -^Df{x) + ^ f:[/^i"^(x) cos(^) + UP{X) sin(!|)]. 



* n=l 

Fields even under the Z2 symmetry will have zero modes and they correspond to 
the SM particles in our usual 4D world whereas fields odd under Z2 symmetry 
will only have KK modes and will be absent in the low energy spectrum so that, 
in the four dimensional Lagrangian, we get rid of the the zero modes of wrong 
chirality (i.e., Qr, Ul, and Dl) and the fifth component of the gauge field, A5. If 
we could measure and study their properties, they would tell us everything about 
the higher dimensional space. 

The elecroweak, SU{2)l x t/(l)y, part of the SM Lagrangian in 5D can be 
written as follows 

^ = I dy{Cl^^,^^ + £^ + £Let.c + ^''), (4.21) 

where y = x'^ is the coordinate along the extra dimension. The fermionic piece 
of the 5D Lagrangian is defined as 

^Leuc = Q{ff''DM)Q + U{iC''Dm)U + D{iC''Dm)D, (4.22) 
where, M,N = 0,1,2,3,4 corresponds to the five-dimensional Lorentz indices. 
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The fermionic fields Q, U and D denotes five dimensional generic quark doublet, 
up type quark singlet, and down type quark singlet, respectively. Unlike in the 
SM, these fermionic fields have both chiralities and are all vector type. The 
covariant derivative. Dm is defined as Dm = Om — W^mT"' — ig'BMY with g is 
the five dimensional gauge coupling constant of the group SU{2)l and g' is that 
of f/(l)y. Here, T'^ and Y are the corresponding generators. In addition, Tm 
are the five dimensional gamma matrices with T ^ = 7^ and r4 = ^75. The Higgs 
piece of the 5D Lagrangian is 

C"" = {Dm^)Kd''<P)-V{<P), (4.23) 

and the gauge piece is 

^kinetic = ~4 Ft^^PuN " -^F^^^Fmn- (4.24) 
j=i 

The field strength tensors associated with the SU{2)l and f/(l)y gauge group 
can be expressed as 

n/TV = duW}, - OnWI, + ge'^^Wi,W^^, (4.25) 

and 

Fmn = OmBn - OnBm, (4.26) 
respectively. Finally the Yukawa piece reads 

Z:^ = gF„</.^f/ + QYd(j)D. (4.27) 

In the above equation, the fields and (f)'^ = ir'^cf)* stand for the standard Higgs 
doublet and its charge conjugated field as each refers to each in order. Finally, 
Yu and Yd correspond the Yukawa matrices in the five dimensional theory which 
are responsible for mixing different generations. For simplicity lepton or gluon 
indices are not included. 

Substituting the KK decompositions of the five dimensional fields given above 
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into the eq. 4.21 and integrating over the extra dimension y, we obtain the effec- 
tive four dimensional Lagrangian. One can reahze simply that the KK excitations 
receive mass not only due to the vacuum expectation value of the zero-mode Higgs 
but also from the kinetic energy term in the five dimensional Lagrangian. The 
mass of the v}^ KK particle is given by = \/mQ + m^. Here mo corresponds 
to the zero mode mass and m„ = n/R. 

Depending on the underlying fundamental theory, extra dimensions may or 
may not be accessible to all fields in the model. According to this accessibility, 
extra dimensions can be grouped into two, including "universal extra dimensions" 
(see for example (UED) [22, 23, 24] and the references therein) and "non-universal 
extra dimensions" (NUED) (see for example [25, 26] and the references therein), 
respectively. In a theory with universal extra dimensions, all fields in the model 
feel the extra dimensions. Conservation of extra dimensional momentum leads to 
the key future of such theories that KK number at each elementary interaction 
vertex is conserved. As a result of this feature, production of an isolated KK 
particle at colliders is forbidden. Instead they are produced in pairs. This, in 
turn, implies that there is no tree-level contribution to weak decays of quarks and 
leptons. They enter into the calculations only through loop corrections. However, 
in a theory with non-universal extra dimensions, some of the SM fields are con- 
fined to a 4D brane and the others live in the bulk. In this case, the Lagrangian 
contains localizing delta function which permits KK number violating couplings. 
Then, the tree-level interactions of KK modes with the ordinary particles can 
exist. 

4.3 The Randal- Sundrum Model 

The large extra dimensions which are discussed in the previous section took the 
advantage of the fact that branes could trap particles and force but neglected 
the energy that the branes themselves could carry. However, according to Ein- 
stein's theory of general relativity, gravitational field is induced by means of 
energy, which means that when branes carry energy, they should curve space and 
time. Lisa Randall and Raman Sundrum [30, 31] tried to explain how space-time 
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would be curved in the presence of two 4D energetic branes that bounded the 
extra dimension of space where the bulk^ geometry is anti-de Sitter^, by solving 
Einstein's gravity equations based on the assumption that both the bulk and and 
the branes have energy. In this space-time the 4D branes and any single slice 
along the fifth dimension are completely flat. But the 5D space-time under con- 
sideration is nonetheless curved. The technical term for this type of geometry is 
'warped'. This section focusses on a warped five- dimensional world that provides 
an alternative approach to explain the huge discrepancy between niEw and Mpi 
without the need for a large extra dimension. In this scenario (called as Randall- 
Sundrum (RSI) Model), the geometry contains two 4D fiat branes, the Planck 
brane where the gravity is localized and the TeV brane where all SM particles are 
confined, that bound a fifth dimension which is compactified to / Z2 orbifold. 
The two 4D fiat branes with opposite tensions, which reside at the orbifold fixed 
points together with a finely tuned non-vanishing 5D cosmological constant A, 
serve as sources for 5D gravity. Since the two branes are completely fiat, the 
induced metric at every point along the extra dimension has to be the ordinary 
fiat 4D Minkowski metric, and the components of the 5D metric depend only on 
the fifth coordinate, y. Thus, the most general space-time metric satisfying these 
properties is given by 

ds"^ = e-^^y'>r]^^dx^'dx'' - dy^, (4.28) 

where e""^^^-* is called as the warp-factor which determines the amount of cur- 
vature along the extra dimension. It is also possible to write this metric in a 
conformally flat form where there is an overall factor. To go into the conformally 
flat frame, all we need to do is to make a coordinate transformation of the form 
z = z{y) such that dy and dz are related by 

g-A{.)/2^^ = dy. (4.29) 

^FuU highcr-dimcnsional space. 

■^Space-time with constant negative curvature. 
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Then, the metric in the eq. 4.28 becomes 

ds^ = e-^^'\7]^^dx^'dx'' - dz^), (4.30) 

and we get 

9MN = e-^^'^-QMN. (4.31) 
where qmn is the fiat metric, qmn = Vmn- 

4.3.1 The Einstein tensor and the brane tensions 

The starting point is the action (see [49] for example) 

S = J d'x^{MlR), (4.32) 
where is the 5D Planck scale, R = g^^ Run and Rmn reads 

RmN = '^MK,N ~ ^MN,K ~^ ^MN^KL + ^ML^NK^ (4.33) 

with the connections 

^MN = g^^^LMN = -g^^{gLM,N + gLN,M — gMN.L)- (4.34) 

After the variation of this action with respect to gMN (see [48] and [49] for details) 
we get 

5S = j d'^xid^MlR + ^Mldg^^'RMNl (4.35) 



where 



= \59MNg''''^g ; ^/'^ = -/^^^^^<5^i.L. (4.36) 



Finally we obtain 

6S = -Ml j d^xiR''"" - \g''''R)^g6gMN. (4.37) 
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with the Einstein tensor Gmn = Rain — \gMNR- Now, we will calculate this 
tensor for the special metric given in eq. 4.31. From now on, since Qmn = Vmn, 
all covariant derivatives Vm (which are with respect to the metric g) will be 
replaced by the normal derivative Om- Let us calculate Rmn term by term. The 
first term is: 

= 2{9^'^^l9sM,K + gsK,Ai — 9mk,s]},n 
= ^{e^V^^[{e-^VsM),K + {e-^VsK),M - (e" V/i^),s]},7v 
= ^{v^'^l-VsAidxA - VsKdAiA + rjAiKdsAlj^N 
= ^{v^'^i-VsAidNdxA - risKdNdAiA + riMKdNdsA]} 
= -^VM^NdxA - ^Vs^NdAiA + ^r]fjdNdsA, 



(4.38) 



where 



Vm = hiN vZ = (4.39) 



with the space-time dimension d. Substituting these equations into Tfjj^ we 
obtain the first term in Rmn as 

^MK,N = —O^OmA. (4.40) 
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By using the same procedure we get the other terms as: 
1 



1 



-^{9 [9sM,N + 9sN,M — 9mn,s]},k 



-{e^V^^[{e-^r]sM),N + {e~^VsN),M - {e~^VMN),s]},K 
^{v^'^l-VsAidNA - VsNdAiA + r]MNdsA]}^K (4.41) 

-^Vm^kOnA - K^^OkOmA + ^v^'^VMNdxdsA 
-dMdNA + ^?7m7v<9^^, 



^^N^KL — 9^^^ SMNQ'"^^ PKL 



KS r 1 LP r 1 

2 fi' [9SM,N + fi'5Af,M — fi'MAf,sJ [fi'Pii',^ + fi'Pi,i^ ~ 9KL,p\ 

^e^Tl^^[{e-'^r]sM),N + (e-Viv),M - {e-^TlMN),s\ 

.LP\_ 



= ^^^'^[-?75m57V^ - VSnQmA + ?7MAr(95A] 

X ?7^^[-?7Pi^9LA - riPLdxA + ?7xl<9p^] 
1 



(4.42) 



-[-ritidNA - t^'^SmA + Ti'^'^TiMNdsA] 
= '^duAdNA - '^7]MN{dA)\ 
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^ML^NK — 9 ^SMLQ ^PNK 

1 r 1 1 J^p r 1 

= 2 5' [9sM,L + gsL,M — 9ml,s\29 [9pn,k + 9pk,n — 9nk,p\ 

= ^ A^^[(e"VAf),L + (e-Vi),*/ - (e"VfL),s] 
X eV^[(e"^r/P7v),it + {e-^r]pK),N - (e-^T/7VK),p] 

X r]^^[-r]PNdKA - VpkQnA - r]NKdpA] 
= ^I-VmOlA - tj^OmA + rj'^'riMLdsA] 

A L A . LP^ 



X [-VnOkA - VkOnA + t'^iqNKdpA] 
= \[-VMdLA - v^OmA + VMLd^'A] 
X [-v^dKA - ti^OnA + tinkOlA] 
= \{2 + d)dMAdNA - ]^T^MN{dAf. 



Finally Rmn becomes 



Rmn = "^-^OmOnA + "^-^OmAOnA 

d — 2 1 
H —VMNidAf - -rjMNd'^A. 



Now, we will find R using R = g Rmn as follows: 

R = c^v^'^^Rmn 

= e^l^LZ^d'A + + - p^A] 

= e^[(l - d)d'A + (2-^Kl-^) (g^)2]^ 

Substituting the eqs. 4.44 and 4.45 into the Einstein tensor we get 

2 — d 1 

Gain = -^—{-dAjAdNA + OmOnA 
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Here d = 5 stands for the number of dimensions. Using this expression we can 
evaluate the non- vanishing terms 6*55 and G^jy. Let us start with G55: 

G55 = —{Ia" + A" - rj,,[-A" + i/l'^]}, (4.47) 

where A' = d^A and 7755 = — 1 and we get 

G55 = —A''. (4.48) 

Gfj,iy can be obtained in the same way: 

G^. = -^{\d^Ad,A + d,d,A - Ti^^id^'dKA - ]^d''AdKA]}. (4.49) 

Since A = A{z), the terms including 4D differentiation vanish and we obtain, 

G,. = lr^,A-A"+^-A''). (4.50) 

Now, we consider the 5D Einstein action for gravity with a bulk cosmological 
constant A: 

S = j Sx^{MlR-K). (4.51) 
Taking the variation of the action with respect to metric 

5S = j d'x[Sy^M^R + ^Mldg^'^^RMN - ^v^A], (4.52) 
we get (see eq. 4.37 for the variation of the first term in eq. 4.51) 

5S = j d'x[-M!G^''' - Ug^''']^5gMN. (4.53) 
Then, we can simply write the Einstein tensor as 

Gmn = -^^9mn- (4.54) 
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The 55 component of Einstein tensor will then be: 



= 4:rA^?55 = \T^e-^^'\ (4.55) 



and y4' reads, 



There exists a solution if and only if A < 0. This means that the important case 
for us will be considering anti-de Sitter spaces, that is the spaces with negative 
cosmological constant. Now, let us take / = e"^*-^-'''^ and, therefore, /' = \A'{z)f . 
Substituting them into the eq. 4.56 we obtain 



/' 1 -A 



(4.57) 



P 2\j 3Mf 

Solving this differential equation we get a relation for / such that, 

/ = e-^(^) = ] (4.58) 
{kz + coY 

where k"^ = If we choose e'^^^'^ = 1, 

e-^(°) = ^ = 1, (4.59) 

^0 

we get Co = 1 and, finally, we obtain 

This solution must be symmetric under z — —z reflection since we are on /Z2 
orbifold and, therefore, we take 
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The RS metric is then obtained as 



ds^ = — — ^ — -iri^ydx^dx" - dz'^). (4.62) 

ik\z\ + ir 



Using the eq. 4.29, we obtain A{z) as follows: 

dz 



k\z\ + 1 

and solving for y we get 



dy, (4.63) 



ln(A; 2; 


+ 1) 


kz/\z\ 





y = , / I + ci. (4.64) 



If we choose y = to correspond to 2; = the constant Ci becomes Ci = 0. Then, 
we have 

^ e-^'^l^l, (4.65) 



{k\z\ + iy 

and the final form of the RS metric in y coordinates becomes 

ds^ = e-^^^y^ri^^dx^'dx" - dy^, (4.66) 

where y corresponds to the physical distance along the extra dimension, since in 
that metric there is no warp factor in front of dy^ term and the Planck (TeV) 
brane is located at y = {y = tq). 

At this stage, we would like to check whether the AD components of the 
Einstein tensor given in eq. 4.50 are also satisfied or not. Using the eq. 4.61, 
A{z) is obtained as 

A{z) = ln{k\z\ + If. (4.67) 
The derivative of A{z) with respect to z gives 



where = e{z) = {0{z) — 6(—z)) and one more derivative of A{z) with respect 
to z reads 
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Substituting A' and A" into eq. 4.50 we get 



G 



3 f 4P 



Ak[5{z) - 5{z - zi)] 
k\z\ + 1 



(4.70) 



Here the first term is the contribution of the bulk cosmological constant into the 
energy momentum tensor. The remaining delta functions should be compensated 
by the additional sources onto the branes. To do this we need to find the enenrgy- 
momentum tensor of a brane tension term V using the action 



S = j d^xVy/^ = j d^xV^=6{y) 



(4.71) 



for a fiat brane at y = 0. Taking the variation of this action with respect the 
metric gfj,y we get 

6S = J d^xV^g^Jo'^'-y^,, (4.72) 
and by using the energy-momentum tensor for a brane 

1 6S 



T, 



we obtain T^^, as 



1 



- fJ,U 



-r9t,vV8{y). 



(4.73) 



(4.74) 



Therefore, at the branes we need to have two brane tensions to compensate the 
delta functions. Thus, we need the equality 



Ak{5{z) - 5{z - zi)) 
k\z\ + 1 



Vp^v \ Vo5{z) + V^5{z- zi) 
2M3 [ k\z\ + 1 



(4.75) 



Then, one can simply conclude that, the brane tensions at the two fixed points 
will have opposite signs and they are given by 



Vo = -V, = UkMl 



(4.76) 



Finally, using the expression k"^ = the bulk cosmological constant is obtained 
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as 



A 




(4.77) 



12M3 



4.3.2 The Radion 



There is a potential gap in this scenario that needs to be filled. In the Randall- 
Sundrum scenario, it is assumed that the brane dynamics would naturally lead to 
branes to be located at a modest distance apart but it is not explicitly explained 
that how the distance between the two branes is established since their solution 
is valid for any choice of tq. If the distance between the two branes remains 
undetermined, when the energy or the temperature of the universe evolves, the 
branes will have the potential to move toward or against to each other. If the 
brane separation could change, the universe would not evolve in the way it is 
supposed to in 4D and thus the warped 5D universe would not agree with the 
cosmological observations. Goldberger and Wise (GW) [50] did the important re- 
search that closed this gap in the theory by fixing^ tq ~ 30/A; without introducing 
any large finetuning. They suggested that, in addition to the graviton, there is a 
massive particle that lives in the 5D bulk for which the equilibrium configuration 
for the field and the branes would involve a modest brane separation. Denoting 
the scalar field in the bulk by 0, the action under consideration will be (see for 
example [47] and [48]) 



where the first term is the usual 5D Einstein-Hilbert action and the second term 



the induced metric on the branes, are the brane induced potentials for the scalar 

field on the Planck and on the TeV branes. We will look for an ansatz of the 

^In the negative tension brane where we Hve all mass scales are exponentially suppressed, 
e-^''«Mpi - 1 TeV. Therefore, kr^ ^ InlQ^^ - 30 since the Planck scale Mpi = 10^^ TeV. 




(4.78) 




is the bulk action for the scalar field, while the next two terms, with being 
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background metric again of generic form as in the RS case such that, 



_ e-^^^y^r]^^dx''dx'' - dy^. (4.79) 

The Einstein equations will be exactly the same as we have derived for the RS 
case, except the energy momentum tensor that is derived from the action of the 
scalar field. Let us take the variation of this action with respect to metric term 
by term. The variation of the first term is given in eq. 4.37. The metric variation 
for the second term 

rf'xv/^[^V0V</. - vm = j rf^a;v/^[^^7^^^9M09^0 - (4.80) 

reads 

5 j d^a:v/^([^V0V0- V(0)] = j d'>x5^[\{dct>f - V(0)] 

d''x^\5g^'''dM<i>dM4> 

d'x{-^g^HgKL)[-{d4>f - vm 
<fx^]^{g''''g''HgKL)dM(^dr,(^. 

By making simplification, we obtain the variation as, 

5 j d'^x^iU/ct^Vct^-Vm = J d'x{^-g^'^[^-{d<py-Vm 

-^-d^^d'^^j^SgKL. 



(4.82) 



Finally, we will look at the metric variation of the last two terms which can be 
written in a more compact form as 

j rf'xV^Ap(0) + j rf^xy^AT(0) = j d'x^^ Yl - ^0, (4.83) 
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where i denote TeV and Planck branes. If we do the metric variation to the 
combination of third and the fourth terms we get 

j d''x^,Xp{4>) + J d^'x^^Xri^) = J c^'x^v^-^ J] A,(</.)5(2/ - y,) 

J ^ Vfi'55 , 



i 

(4.84) 

Replacing K and L with M and N, respectively, we obtain T^^^ as 

1 1 ^ (4-85) 



By considering the equality (see Appendix B for its derivation) 

Rmn = f^'^TMN, (4.86) 

where 

Tmn = Tmn — -gAmT, (4.87) 

we get 

4^'2 _ ^ -^V{<Po) - y 5^ K{<Po)5{y - Vi). (4.88) 

i 

On the other hand, the 55 component of the Einstein equation is obtained as 

G'55 = y0o'-«V(0o), (4.89) 



and we get 



A" = ^<l>'o'-jV{<Po). (4.90) 
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In these equations (po denotes the solution of the scalar field, which is assumed to 
be only a function of y: (p = (poiv)- In addition to these two equations, the bulk 
scalar equation of motion is found by using 

^^7^^ - l§ = 0' (4.91) 

d{dM(j)) d(j) 

as 

f)V 

5a/(v^/'^9^0) = -g^V9- (4.92) 
We can write the term dM{\/gg^^^ dN(t>) also in the form 

(4.93) 

Since (p = 4>o{y) and A = A{y) we get 

fiv 

-e-^^<j)'^ + 4e-^^A>'o = -^e-^^. (4.94) 



Including the brane tensions we can write 



,// . ./ ,/ dV{(po) ^ dXi{(f)o) 

6o-4A0o = — — + 6{y-yi). (4.95) 



The metric itself have to be continuous. However, there is no requirement that 
the derivative of the metric to be continuous. In the first Einstein equation given 
in 4.88 there exists the explicit delta function term —^"^i-pT K{<Po)S{y — yi) 
at the branes. It seems not the be balanced by anything else unless there is a 
jump in the derivative A' at the branes. If the derivative jumps from A'(0 — e) 
to ^'(0 + e), this implies that locally A' contains a term of the form 

^[y = 0) ~ [A'(0 + e) - A'(0 - eMy), (4.96) 
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where e{y) is the unitstep function. Taking one more derivative with respect to 
y we get 

A"{y = 0) ~ [A'iO + 6) - A'{0 - e)]S{y). (4.97) 

Thus, the delta function is proportional to the jump of the derivative of A'. In 
the same way, the delta function in the bulk scalar equation of motion will be 
proportional to the jump of the derivative (pQ. Therefore, the boundary conditions, 
(or jump equations) will be given by 

[A']. = 



These are coupled second order differential equations. So, they are quite hard 
to solve. Only for specific potentials solution can be simplified. Defining the 
function W^cj)) such that 



A' = -W{<Po), 



and substituting these equations into the 55 component of Einstein equation we 
get 

n0) = ^(^)^-yH^(0f. (4.100) 

This is called as the consistency equation since when we plug in the expressions 
for A' and 0q into the Einstein and scalar equations, we will find that all equations 
are satisfied. Then, the jump equations are given by 

1 dW _ dX,{<Po) ^^-^^^^ 

2 ^ 06 ^' 86 
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If W were given, the coupled second order differential equations would be reduced 
to first order equations that are easy to solve. We can simply start with a super- 
potential that will produce V with the required properties. In our case, we would 
like the bulk potential to include cosmological constant term (independent of (p) 
and mass term (quadratic in 0) in its simplest form. So we choose 



6k 



W{(f)) = — - u(j)'. (4.102) 

where the first term is just one needs for cosmological constant and the second 
one is for the mass term. Then, 

1 dW 

^' = ^a^ = -*. (4^103) 

Substituting (j) = Ce"^^ we get m = —u. If we use the boundary condition that 
at y = 0, = 0p, we find that C = (pp. Then we have 

(f) = (Ppe-'^y. (4.104) 
From this value of the scalar field at the TeV brane at r it is determined to be 

0^ = ^pe""''. (4.105) 
This means that the radius is no longer arbitrary but given by, 

r = -In—. (4.106) 

This is the GW mechanism. The background metric will then be obtained from 
the equation 

A' = ^W{<Po) = k - ^<Ple-'-y, (4.107) 

given by the solution 

A{y) = ky + ^e-'^y, (4.108) 
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where the first term is the usual RS warp factor, while the second term is the 
backreaction of the metric to the non-vanishing scalar field in the bulk. Since 
we want to generate the right hierarchy between the Planck and Weak scales we 
need to ensure that kr ~ 30 from which we get 

^ln(^)~ 30, (4.109) 

which is the ratio that will set the hierarchy in the RS model. Since 4'p/4>t is 
constant, so u is kept constant. 

4.3.3 The coupled field equations and the radion mass 

Once we have established the mechanism for the stabilization of radion, we realize 
that the radion is no longer massless. Then, the question is what will be the value 
of the radion mass in the GW stabilized RS model? For this we consider spin-0 
fluctuations of the coupled gravity-scalar system over the background. This can 
be parameterized in the following way 

ds^ = e-2^(^)-2^("'^)r/^,rfx'^rfx'^ - (1 + ^^^^ 
(j){x,y) = 0o(y) + <^(x,i/). 

It looks like as if there would be three different fluctuations, namely F, G, (f. We 
will use this ansatz to linearize Einstein and scalar field equations. Then, some 
coupled equations for F, G, and ip will be obtained. The linearized Einstein eq. 
B-4 (see Appendix B for details) equations read 

6Rmn = ^^^STmn, (4.111) 
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and by using this equation together with the ansatz given in eq. 4.110 where 



a - p-'2A(y)-2F(x,y) 
au ^ ^2A{y)+2F{x,y) nu 

(4.112) 

(755 = -{1 + G{x,y)f, 
g'' = -a + Gix,y)r', 

the hnearized form of the Rfj_u is obtained as (see Appendix B for the detailed 
calculations) 

SR^, = -2d^d,F + d^d.G - ri^^ap + e-^^r]^, 

X [F" - 8A'F' - AG' - 2FA" + 8FA'^ - 2GA" + SGA'^]. (4.113) 

Inspecting the equation, we realize the dfj_di, term must vanish since to linear 
order in perturbations all the terms in T^^, (see eq. B-61) are proportional to rj^i,. 
Then, one can immediately conclude that G = 2F since 

5R^, = - 2d^d,F + d^d.G + (4.114) 

Substituting this into the equation for 5/2^,^ we get, 

SR^^ = -7]^^UF + r]^^e-'^^{F" - lOA'F' - + 2AA''^F]. (4.115) 

Similarly, 5R^^ and 5R^^ is obtained as 

5i2^5 = -39^F' + QA'd^F, (4.116) 

and 

(5i?55 = -AF" - 2e^^UF + IQAF'. (4.117) 



55 



(see Appendix B for details). Now, we will present the linearized source terms 
5T^y, 6Tfj^5, and 5T55 (see Appendix B for their derivations): 

Sf^u= ^e-2Vb^'(0o)-2^(0o)i^] 



i 

1 



6T,,= --<P'^d,ip, (4.118) 
5T55= -^y(0o)F-i^l^'(</)o)-^>'o 

-|E[^^ + 2A.(0o)F]%-y.). 



Finally, we present the linearized scalar field equation for completeness: 

dcj) 



e Dy-^ +4A^ +— -^^ = -}_^( y + — F)%-yO 



(4.119) 



Using the equation SRn5 = i^'^ST^^ we get 



3{d^F' - 2A'd^F) = K^(P'od^ip. (4.120) 



This can be integrated immediately to obtain 



hF'-2A'F) + C{y), (4.121) 



k2 



where C{y) is the integration constant. We set this constant as zero since we 
require that the fluctuations F and (p are also localized in x. Let us find the 
boundary conditions for F and if on the two branes. Now, we use the equation 
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= f^'^^T^b which leads to, 



and, taking in to account the continuity of the metric but not its derivative, we 
get the equation 

ri = f^A.(^o)F + ^^V,. (4.123) 

Here the delta function is proportional to the jump of the derivative of F', 

F"{y = 0) ~ [F(0 + e) - F'(0 - e)]5(2/), (4.124) 

with 

F\y = 0) ~ [F'(0 + e) - F'(0 - e)\e{y), (4.125) 

where £(?/) is a unitstep function. In the same manner, by using linearized scalar 
field equation, eq. 4.119, we can take 

Using the jump equations, eqs. 4.98 and 4.121, with C{y) = we get 

= — M0o)F + -^^v^. 

Thus, it provides no new constraints. Then, only the second boundary condition 
must be taken into account. For a convenient limit ^ ^'^2'°'' ^ 1? the second 
boundary condition is simply = 0. Then, in this limit the first boundary 
condition is just 

{F' -2A'F)\i = 0. (4.128) 



57 



Now a single equation for F is obtained as follows. Considering the combination 
e^^SR^u + Tj^ydR^^ in the bulk we get 

e^^5R^, + r/^.5i?55 = Sr/^.fe^^DF + F" - 2AF' + 2F{A' - (4.129) 

where A' — 4^4'^ = '^V{(j)Q). Here, we do not take into account the 5 terms since 
we work in the bulk. The similar combination for the source terms reads 

e'^K'df^, + r^^^/c^^fgs = Ak\,,V{<P^)F + 2k'7]^,(I)'^^' , (4.130) 

and equating them to get 

e^^DF + F" - 2A'F' = — </.'o<^'. (4.131) 

3 

Using eq. 4.121, (f' can be obtained in terms of F as 

3 F' - 2A'F 

^ = <^-^=''> 

If we take one more derivative with respect to y we get 

3 [F" - 2A'F' - 2A"FW, - {F' - 2A'F)<j,'^ 
^ = ' ^^-'^^^ 

and by making some arrangements in this equation we obtain, 

= 2{F" - 2A'F' - 2A"F) - 2{F' - 2A'F)^. (4.134) 

3 00 

If we substitute into the eq. 4.131 we get 

F" - 2A'F' - 4:A"F - 2^F' + AA'^ = e^^DF, (4.135) 

to be solved in the bulk. It is important to note that each eigenmode DF^ = 
—m'^Fn to this equation has two integration constants and one mass eigenvalue. 
The first constant corresponds to the overall normalization, while the remaining 
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one is fixed by the boundary condition on the Planck brane, and the mass is 
determined by the boundary condition on the TeV brane. 

Now we are ready to calculate the radion mass. In the following we will show 
how backreaction generates a non-vanishing mass for the radion field by using eq. 
4.135. Substituting 0o(y) = 0pe~"^ into the above equation we get 

F" - 2A'F' - AA"F + 2uF' - AuA'F + m^e^'^F = 0, (4.136) 

where A(y) is given in eq. 4.108. The backreaction will be treated as perturbation 
such that 

A{y) = k\y\ + ^^e-'-^y\ (4.137) 



where / = /t0o/v^- Then, we will look the solution in terms of the perturbative 
series in /. The solution is written as follows 

Fo = e''\y\{l + Pfo{y)) ; ml = I'm'. (4.138) 

By substituting the solution into the eq. 4.136 and keeping only the leading terms 
in P we get 

+ 2{k + m)/^ = -mV'lJ^I - - M)Me-2"l^l. (4.139) 
By solving this equation we get 

M = Ce-'^'^^-^\y\ - ,^0—^e'^\y\ - 2^-^^e-'^\y\. (4.140) 

If we make the same substitution in the boundary condition F' — 2A'F = we 
get 

/2/^ + lufe-'^'^y^ + ^uZVoe-'"l^l = 0, (4.141) 
and by keeping only the leading terms in 1' we obtain 

/o + \ue-'''\y\ = 0. (4.142) 

at the boundary, where \y\ = tq. If we substitute this into the eq. 4.139, m is 
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obtained as 



= 1^^^2g-2(.+fc).o^ (4143) 

3 ru 



and the radion mass reads 

mlad^on = 3^ e ^ °- (4.144) 

4.3.4 Coupling to SM fields and the normalized radion 
field 

In the previous section, including the backreaction, a mass scale 0{TeV^) is 
obtained for the radion mass. Then, the wavefunction can be written as 

Fo(a;,i/) = e2%l(l + /2/o(2/)W, (4.145) 

where fo{y) obtained using the integral in the eq. 4.140. Based on the as- 
sumption P <^ 1, we see that the backreaction induces a small correction to 
the unperturbed wavefunction. So for purposes of determining the coupling of 
the radion to the TeV brane, we can use the unperturbed wavefunction, namely 
F(x,y) = e^'^^y^ R{x) . Let us try to find the coefficient of {dF^ term in the La- 
grangian ^/gR with ^/g = e~^^~'^^{l + 2F), so that we are able to write the 
normalized wave function R{x). R can be splited as i? = g^'^R^p + g^^R^^. Then, 
using the R^y (see Appendix B for details) read 

d d F -2A-2F 

R,, = -2d,d,F + 2-^^-r],,nF + r]^, 



X [A" + F" - 4{A' + F 



1 + 2F ' '^"(1 + 2^)2 

2{A' + F')F' 



2F) 



^tt¥" " ^^'^^rri^ " ^^'-^^'^^ + ^vAm'- (4.146) 
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and multiplying by ^/gg^^ from left, for the first term in R, including [dF)'^, we 
get 

^g'^'R^^ ^ e~2^(l - 2F)(-4(1 + 2F)FUF - 4(1 + 2F)UF + 6(1 + 2F){dFf 
-^dFf...), 

(4.147) 

where FBF = Fd^d^'F = d^,{Fd^F) - {dFf. Substituting this quality into the 
above equation we obtain, 

g^'^R^^ ^ e-2^(l - 2F){-4[d^{Fd^F) - {dFY]{l + 2F) 

- A{d^{Fd^'F) - {dFf)F - 8[9^,(F9^F) - [dFf] (4.148) 
-2{dFf + ...). 

Similarly, using the equation for R^^ (see Appendix B for details) below 

i?55 = -4(A" + F") + 46^^+2^(1 + 2F){dFf 

A' + F'] 
1 + 2F 



A' J- P'^ P' (4.149) 
2^2A+2F^p^^ + 2F) + ^ J + AiA' + F'f, 



the second term in R is obtained as 

\/gg ^55 (1 + 2F)2 

_ ^ ^ 4(A- + n _ 4e^-^+^^(gF)^ 2e^^+^^DF 

^(1 + 2F)2 1 + 2F + 1 + 2F 

8(A' + F')F' ^ 4(^' + F')^ 



:i + 2F)3 (1 + 2F)2 

(1 - 2F) [a{A" + F"){1 - AF) - 4e2^(l - AF^){dFf 

26^^(1 - AF^){d^,{Fd^'F) - {dFf)F - 8(A' + F')F\l - QF) 
4(A' + F')'(1 -4F)). 

(4.150) 
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So, one can simply find the coefficient of [dFY in y^R in linear order as e ^"^6. 
Then, a straightforward calculation gives 

-Ml j dy^R D 6Ml{dRY J e'^^e^'^l^l = ^(e^'^^'^" - l){dR)''. (4.151) 

Therefore, the normalized radion r{x) is R{x) = r{x)e~'''^^ / \/6Mpi, which is 
obtained by using M^/k = M^J2. 
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CHAPTER 5 



LEPTON FLAVOR VIOLATING 
RADION DECAYS IN THE 
RANDALL-SUNDRUM 
SCENARIO 

The hierarchy problem between weak and Planck scales could be explained by 
introducing the extra dimensions. One of the possibility is to pull down the 
Planck scale to TeV range by considering the compactified extra dimensions of 
large size [20, 21]. The assumption that the extra dimensions are at the order 
of submilimeter distance, for two extra dimensions, the hierarchy problem in the 
fundamental scales could be solved and the true scale of quantum gravity would 
be no more the Planck scale but it is of the order of EW scale. This is the 
case that the gravity spreads over all the volume including the extra dimensions, 
however, the matter fields are restricted in four dimensions, so called 4D brane. 
Another possibility, which is based on the non-factorizable geometry, is introduced 
by Randall and Sundrum [30, 31] and, in this scenario, the extra dimension is 
compactified to /Z2 orbifold with two 4D brane boundaries. Here, the gravity 
is localized in one of the boundary, so called the Planck brane, which is away from 
another boundary, the TeV brane where we live. The size of extra dimension is 
related to the vacuum expectation of a scalar field and its fluctuation over the 
expectation value is called the radion field (see section 2 for details). The radion 
in the RSI model has been studied in several works in the literature [50]- [58] (see 
[47] for extensive discussion). 

In the present work, [63] we study the possible LEV decays of the radion field 
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r. The LFV interactions exist at least in one loop level in the extended SM, so 
called uSM, which is constructed by taking neutrinos massive and by permitting 
the lepton mixing mechanism [41]. Their negligibly small Brs stimulate one to 
go beyond and they are worthwhile to examine since they open a window to test 
new models and to ensure considerable information about the restrictions of the 
free parameters, with the help of the possible accurate measurements. The LFV 
interactions are carried by the FCNCs and in the SM with extended Higgs sector 
(the multi Higgs doublet model) they can exist at tree level. Among multi Higgs 
doublet models, the 2HDM is a candidate for the lepton flavor violation. In this 
model, the lepton flavor violation is driven by the new scalar Higgs bosons S, 
scalar and pseudo scalar A^, and it is controlled by the Yukawa couplings 
appearing in lepton-lepton-S vertices. 

Here, we predict the BRs of the LFV r decays in the 2HDM, in the framework 
of the RSI scenario. We observe that the BRs of the processes we study are at 
most of the order of 10~^, for the small values of radion mass and their 
sensitivities to decrease with the increasing values of m^. Among the LFV 
decays we study, the r ^ /i^ decay would be the most suitable one to measure 
its BR. 

5.1 The LFV RSI model radion decay in the 
2HDM 

The RSI model is an interesting candidate in order to explain the well known 
hierarchy problem. It is formulated as two 4D surfaces (branes) in 5D world in 
which the extra dimension is compactified into S^/Z2 orbifold. In this model, 
the SM fields are assumed to live on one of the brane, so called the TeV brane. 
On the other hand, the gravity peaks near the other brane, so called the Planck 
brane and extends into the bulk with varying strength. Here, 5D cosmological 
constant is non vanishing and both branes have equal and opposite tensions so 
that the low energy effective theory has fiat 4D spacetime. The metric of such 
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5D world reads 



ds^ = e-^^^^^r]f,^dx^dx'' -dy\ (5.1) 

where A{y) = k\y\, k is the bulk curvature constant, y is the extra dimension 
parametrized a.s y = R6. The exponential factor e"'^^ with L = Rir, is the warp 
factor which causes that all the mass terms are rescaled in the TeV brane. With 
a rough estimate L ~ 30/A;, all mass terms are brought down to the TeV scale. 
The size L of extra dimension is related to the vacuum expectation of the field 
L{x) and its fluctuation over the expectation value is called the radion field r. In 
order to avoid the violation of equivalence principle, L{x) should acquire a mass 
and, to stabilize r, a mechanism was proposed by Goldberger and Wise [50], by 
introducing a potential for L{x). Finally the metric in 5D is defined as [51]. 

ds^ = e-^My)-2F{x) - (1 + 2 F{x)) dy^ , (5.2) 

where the radial fluctuations are carried by the scalar field F{x), 

F{x) = ^^^^—^r{x). (5.3) 



VQMpi e- 



Here the field r{x) is the normalized radion field (see [48]). At the orbifold point 
9 = IT (TeV brane) the induced metric reads, 

^.n.^^-2^(L)-2..W^^^_ (5.4) 

Here the parameter 7 reads 7 = -^^^ with A = Mpie~^^ and v is the vacuum 
expectation value of the SM Higgs boson. The radion is the additional degree of 
freedom of the 4D effective theory and we study the possible LFV decays of this 
field. 

The FCNCs at tree level can exist in the 2HDM and they induce the FV 
interactions with large BRs. The FV r decays, r lilt^ '^^^ exist at least in 
one loop level in the framework of the 2HDM. The part of action which carries 
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the interaction, responsible for the LFV processes reads 



Y 



(5.5) 



where L and R denote chiral projections L{R) = 1/2(1 =1= 75), 0j for i = 1, 2, are 
two scalar doublets, Ul {Ejn) are lepton doublets (singlets), ^fj ^ and rjfj, with 
family indices are the Yukawa couplings and ^fj induce the FV interactions in 
the leptonic sector. Here (7*"'^ is the determinant of the induced metric on the TeV 
brane where the 2HDM particles live. Here, we assume that the Higgs doublet 
01 has a non-zero vacuum expectation value to ensure the ordinary masses of 
the gauge fields and the fermions, however, the second doublet has no vacuum 
expectation value, namely, we choose the doublets (pi and 02 and their vacuum 
expectation values as 



1 
72 








(5.6) 



and 



< 01 > = 




; < 02 >= 



(5.7) 



This choice ensures that the mixing between neutral scalar Higgs bosons is switched 
off and it would be possible to separate the particle spectrum so that the SM par- 
ticles are collected in the first doublet and the new particles in the second one 
^. The action in eq. (5.5) is responsible for the tree level S — li — I2 {h and I2 
are different flavors of charged leptons, 5* denotes the neutral new Higgs boson, 
S = A°) interaction (see Fig. 5. 3-d, e) and the four point r — S — li — I2 
interaction (see Fig. 5.3-c) where r is the radion field. The latter interaction is 
coming from the determinant factor ^/—g^ = e~'^^^^^~'^'^^^^\ Notice that the 
term e"'^'^'^^^ in \/— is embedded into the redefinitions of the fields on the 



^In the following, we replace with where "N" denotes the word "neutral". 
^Here Hi {H2) is the well known mass eigenstate hP [A^). 
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TeV brane, namely, they are warped as S ^ ^ML) S^arp, I — > e^^^^^^/^ Iwarp and in 
the following we use warped fields without the warp index. 
On the other hand, the part of new scalar action 

^2 = J d^x^/^^lr''^^''{D^(j)2yD,(P2-ml4cf^2] (5.8) 



leads to 



- e-''i''imioh'h' + m'^oA'A^)^, (5.9) 

which carries the S — S — r interaction^ (see Fig. 5.3 1-b). 

Finally, the interaction of leptons with the radion field is carried by the action 
(see [52]) 



J dW-g'""" (^g'""^ hi. iDj-mi n j , (5.11) 



where 



D^l = d,l + ]^wfT.ad: (5.12) 

with Safe = 7 [7a, 7b]- Here w"^ is the spin connection and, by using the vierbein 



■^In general, there is no symmetry which forbids the curvature-scalar interaction, 

S(^^ j d^x^/^g^S^UH^ H , (5.10) 

where ^ is a restricted positive parameter and H is the Higgs scalar field [47, 48, 59]. This inter- 
action results in the radion- (SM or new) Higgs mixing which can bring a sizeable contribution 
to the physical quantities studied. Here, we assume that there is no mixing between first and 
second doublet and only the first Higgs doublet has vacuum expectation value. Therefore, we 
choose that there exists a mixing between the radion and the SM Higgs field, but not between 
the radion and the new Higgs fields. This is the case that the lepton flavor violation is not 
affected by the mixing since the SM Higgs field is not responsible for the FCNC current at tree 
level. 
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fields e^, it can be calculated (linear in r) as 

w;' = -^d,r{e'^'e';-e^^e''^), (5.13) 

(see Appendix C for details). Notice that the vierbein fields are the square root 
of the metric and they satisfy the relation 

Using eqs. (5.11)-(5.14), one gets the part of the action which describes the tree 
level I — I — r interaction (see Fig. 5.3-a) as 

5'3 = fd^x\-?>-rTim-?>—li3rl + A-mirll\. (5.15) 
J \ V 2v V \ 

Now, we are ready to calculate the matrix element for the LFV radion decay. 
The decay of the radion to leptons with different flavors exits at least in one loop 
order, with the help of internal new Higgs bosons S = h^, A^. The possible vertex 
and self energy diagrams are presented in Fig. 5.2. After addition of all these 
diagrams, the divergences which occur in the loop integrals are eliminated and 
the matrix element square for this decay is obtained as 



\M\' = 2[mi - {m,-+m,^y) \A\\ (5.16) 

where 

A fself . pself I rvert , rvert , rvert , rvert /'c; 1'7^ 

^ ~ ho + J AO + JhO + J AO + JhOhO + J AO AO ; [^■^' ) 
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and their explicit expressions are given by 



fself 
JhO 



f 



■self 
AO 



fve 
JhO 



■vert 



7 



128t;7r2 (w'l^-Wh) Jo 

(o I K \ 1 ^hhO^hO 

X i6Wj-—bWhj In 



dx niffi I [iqj {x -l)wh- vt ^ 



ih 



jself 2 



7 



128 V 7r2 {w'j^ - wa) 
X ^3 w'a — 5 wa^ In 



jself 2 
1 AO ''''AO 



{ri( {x-l)w'j^ + r]f ZiAj {hw\-?, In 



jself 2 
^2,A0 '"'AO 



7 



128 V 7r2 



pl pl—x 








TUhO 
T ver 









"if 3 zl^ (y (1 -y)w'f^ + x^{Ay- 1) Wh 



+ x{{l-3y)wh + y{4y- 3) w'^)^ + 5 z^ ((2 y - 1) w', + {2 x - 1) 

- 3{x + y-l)(^xiAx-3)wl + yiAy-3)w'^^ 

- 3whw'^{x + y -1) (^{1 - y + X {Ay - 2)) Wh + {1 - 2y + X {Ay - 1)) w'f, 
+ 3{x + y-l)(^{2x-l)wh + {2y-l)w'^^ 

+ Vt Zih + y - 1) ( - 4 + 2 + w'^ (81/ - 3) + wl (8x - 3) ) 

4 + 2:^h((8 2/-3)a;-3y)) 



' ver ^2 



ni,fi In^^-^ I Qr^f (<(2y-l)+M;42x-l) 
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rvert 
JAO 



fve 



vert 



7 



128 WTT^ 



/•I fl — X 






dx dy < 


T ver 











vl 3 (y (1 -y)w'A + x"^ {Ay - 1) wa 



+ x{{l-3y)wA + y{Ay-S)w'^)') +5z^ (^{2y - l)w'^ + {2x - 1)wa 

- 3{x + y-l) (x{4:X-3)w\ + y{4:y-3)w'^^ 

- 3wAw'A{x + y-l) ({1 - y + X {Ay - 2)) WA + {1 - 2y + X {Ay - 1)) w'a 



+ 3{x + y-l) (^{2x-l)wA + {2y-l)w'^^^ 

+ vfziA i^x + y-l) {^-A + 2w'AWA + wl{%y-3) + w\{%x-3)^ 
- (8z5i + ^'A((82/-3)x-3y)) j 



— ntAo In 



7 



Tver ^2 



dx 



9r]Y(w'A {2y - 1) + wa{2x -1)^ +87]+ ZiA 

vY l4h{y-'^ + x{i - 4y)) {xwh + yw'h) 



dy 



T ver 



64f vr^ 

+ y{x + y-l)w'J{Ax-l)wl + {Ay-l)w'^ 



+ wlx{x + y- 1) {Ax -1) + {x + y - 1) (^2 y ty;, + x u;/, (2 + (4 ?/ - 1)) j 
+ Vt (^i^ + y- 1) ^^ih ((4y - 1) w'^ + {Ax-l)wl + 2] 
- Zihz'^f^(^{Ay-l)x-y + l^^ 

r;^ (wUl - 6 + u;, (1 - 6 x)) - 4 r/+ z,i}j | , 



m^o In ■ 



tier ™2 
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£vert ^ 



fl fl—X 






dx dy < 






r tier 




^ ^AO AO 





+ y{x + y -l)w'j^[{Ax -l)w\ + {Ay -l)w 



+ w^x {x + y - 1) {Ax - 1) + {x + y - 1) (2yw'^ + xwa (2 + w'j (4?/ - 1)) j j 
+ Vt (^{x + y- 1) ((42/ - 1) w'2 + {Ax-l)wl + 2^ 
- ZiAZ^,^(^{Ay-l)x-y + l^^ 

Tver ^2 I , , . I 

_ m^o/n -^° Uf(t^:,(l-6y) + i/;A(l-6x))+4r/+z,^U, (5.1^ 



where 



r se// _i|2'2 |/2 '2 i\ 

-^2,/iO{AO) — J- + a; '"^/i(A) + ^ \^ih{iA) ~ '^h(A) ~ -"-J 



K^{AO) = wl^A) + - 1) {w'h{A) y-l)+x{y w'^i^A) + (y - 1) wl^A) - y ^rh(rA) " 1) > 
^rohO(AO Aof^^ W^h(A) + (1 + W^MA) - 1)) 2/ + X (1 + - 1) + W^i^A) 2/ - 2/) , 

(5.19) 

With the parameters w^^a) = w^^^^ = Zr^irA) = z^hiiA) = 



and 



rrij 



vt — ^N,ili^N,il2 ^ ^N,lii^N,l2i ■ (5.20) 

In eq. (5.20), the flavor changing couphngs i-i represent the effective interaction 
between the internal lepton i, {i = e, /i, r) and the outgoing j = 1 {j = 2) lepton 
(anti lepton). Here, we choose the couplings real. 

Finally, the BR for r ^ /j~ can be obtained by using the matrix element 
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square as 

where is the total decay width of radion r. In our numerical analysis, we 
consider the BR due to the production of sum of charged states, namely 

BR{r If If) = nr^{hh + kh)) ^^22) 

i ^ 

5.2 Numerical Analysis and Discussion 

In four dimensions, the higher dimensional gravity is observed as it has new states 
with spin 2,1 and 0, so called, the graviton, the gravivector, the graviscalar. These 
states interact with the particles in the underlying theory. In the RSI model 
with one extra dimension, the spin gravity particle radion r interacts with the 
particles of the theory (2HDM in our case) on the TeV brane and this interaction 
occurs over the trace of the energy-momentum tensor with the strength l/A^, 

Ani ^ ' (5.23) 

where is at the order of TeV. The radion interacts with gluon [g) pair or 
photon (7) pair in one loop order from the trace anomaly. For the radion mass 
TTir < 150 GeV , the decay width is dominated by r — > gg. For the masses which 
are beyond the WW and ZZ thresholds, the main decay mode is r — WW. In the 
present work, we study the possible LFV decays of the RSI radion in the 2HDM 
and estimate the BRs of these decays for different values of radion masses. We 
take the total decay width F^ of the radion by considering the dominant decays 
r ^ gg (77, //, W~^W~ , ZZ, SS) where S are the neutral Higgs particles (see [55] 
for the explicit expressions of these decay widths). Here, we include the possible 
processes in the F^ according to the mass of the radion. 

The flavor violating r decays r l^lf can exist at least in one loop level, in 
the framework of the 2HDM and the flavor violation is carried by the Yukawa 
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couplings In the version of 2HDM where the FCNC are permitted, these 

couphngs are free parameters which should be restricted by using the present and 
forthcoming experiments. At first, we assume that these couplings are symmetric 
with respect to the flavor indices i and j. Furthermore, we take that the couplings 
which contain at least one r index are dominant and we choose a broad range 
for these couplings, by respecting the upper limit prediction of ^^^^ (see [60] and 
references therein) which is obtained by using the experimental uncertainty, 10~^, 
in the measurement of the muon anomalous magnetic moment and by assuming 
that the new physics effects can not exceed this uncertainty. For the coupling 
^^T-e, the restriction is estimated by using this upper limit and the experimental 
upper bound of BR of /i ^ 67 decay, BR < 1.2 x 10^^^ [61]. Finally, this coupling 
is taken in the range 10"^ — 10^^ GeV (see [62]). For the Yukawa coupling 
we have no explicit restriction region and we use the numerical values which are 
greater than Throughout our calculations we use the input values given in 

Table (5.1). 

Table 5.1: The values of the input parameters used in the numerical calculations. 



Parameter 


Value 


rriho 


0.106 (GeV) 

1.78 (GeV) 

100 (GeV) 

200 (GeV) 
1.1663710-5(Ge1/-2) 



In Fig. 5. 3 we present dependence of the BR (r fi^). The solid- 

dashed lines represent the BR (r ^ r± ^±) for = 100 GeV, = 10 GeV- 
= 10 GeV, ^f ^^^ = I GeV. It is observed that the BR (r ^ r±/i±) is of 
the order of the magnitude of 10"^ for the large values of the couplings and the 
radion mass values ~ 200 GeV^. For the heavy masses of the radion the BR is 
stabilized to the values of the order of 10~^. 

^Thc dimcnsionfuU Yukawa couplings ^§ are defined as ^§ — \ / ^^f- ij- 
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Fig. 5. 4 is devoted to dependence of the BR (r e^) and BR (r — > 

/i^ e^). The sohd-dashed hnes represent the BR (r e^) for = 100 GeV, 

l^^^^ = 0.1 GeV- Ifj^^^ = 10 GeV, ^f^^, = 0.1 GeV. The small dashed line rep- 
resents the BR (r /i±e±) for ^§^^^ = IGeV, = 0.1 GeV. This figure 
shows that the BR (r — ^ r^/i^) is of the order of the magnitude of 10~^^ for 
the large values of the couplings and the radion mass values ~ 200 GeV. For the 
heavy masses of the radion, this BR reaches to the values less than 10~^^. The 
BR (r — > /i^ e^) is of the order of 10^^^ for ~ 200 GeV and for the interme- 
diate values of Yukawa couplings. These BRs, especially BR (r — > /x=''e^), are 
negligibly small. 

Now, we present the Yukawa coupling dependencies of the BRs of the decays 
under consideration, for different radion masses . 

Fig. 5. 5 represents the ^^,tt dependence of the BR (r js^) for ^^.^^ = 

10 GeV. The solid-dashed-small dashed lines represent the BR for the radion 
masses = 200 — 500 — 1000 GeV^. This figure shows that the BR is sensitive 
to the radion mass and, obviously, it is enhanced two orders of magnitude in the 
range 10 GeV < ^§^^^ < 100 GeV. 

In Fig. 5. 6, we present the ^^^^ dependence of the BR (r — * e"*") for ^^^^ = 
0.1 GeV. The solid-dashed-small dashed lines represent the BR for the radion 
masses = 200 — 500 — 1000 GeV. Similar to the r ^ fi^ decay, the BR is 
strongly sensitive to the radion mass. 

As a summary, the LFV decays of the radion in the RSI model strongly depend 
on the radion mass and the Yukawa couplings. The BR for r — > //^ decay is of 
the order of 10~^ for the small values of radion mass and it decreases with the 
increasing values of m^. On the other hand, the BRs for r ^ (r — /i^ e^) 

decays are of the order of 10~^^ (10~^^) for the small values of m^. These results 
show that, among these processes, the LFV r ^ /i^ decay would be the most 
appropriate one to measure its BR. With the possible production of the radion 
(the most probable production is due to the gluon fusion, gg ^ r [55]), hopefully, 
the future experimental results of this decay would be useful in order to test the 
possible signals coming from the extra dimensions and new physics which results 
in flavor violation. 



74 



5.3 The vertices appearing in the present work 

In this section we present the vertices appearing in our calculations. Here S 
denotes the new neutral Higgs bosons and A^. 
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Figure 5.1: The vertices used in the present work. 
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12^ 

(e) (f) 

Figure 5.2: One loop diagrams contribute to r decay due to the neutral 

Higgs bosons Hq and Aq in the 2HDM. i represents the internal lepton, (/^) 
outgoing lepton (anti lepton), internal dashed line the Hq and Aq fields. 
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Figure 5.3: m,. dependence of the BR (r fi^). The sohd-dashed hnes 

represent the BR(r ^ r±/i±) for = 100 Gel^, = 10 GeV- = 

10 GeV, ^1,^ = 1 GeV. 
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Figure 5.5: ^^^-^ dependence of the BR(r — > r^/x^) for = 10 GeV. The 

sohd-dashed-small dashed hnes represent the BR for the radion masses rrir = 
200- 500 - 1000 GeV. 
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Figure 5.6: ^§tt dependence of the BR(r r^'^e^) for ^^^-^ = 0.1 GeV. The 
sohd-dashed-small dashed lines represent the BR for the radion masses rrir = 
200- 500 - 1000 GeV. 
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CHAPTER 6 



CONCLUSION 

Over the last thirty years or so, the SM has been subjected to many diverse 
experiments and most of these experiments have been found to be consistent 
with the predictions of the SM except the recently performed ones on the FCNCs 
including LFV interactions. In the framework of the SM, LFV interactions exist 
in the z/SM by taking the neutrinos massive and permitting the lepton mixing 
mechanism [41, 42] in order to accommodate the present data on neutrino masses 
and mixings. However, the neutrino masses are so small that, the predicted 
BRs of the LFV interactions are too small to explain the experimental data 
obtained. Therefore, it is clear that there must be a theory which goes beyond 
the SM, but which reproduces the results of the SM where the SM has been 
shown experimentally to be correct. The LFV interactions have still been probed 
at many experiments such as the MEG experiment [64], the MEGA experiment 
[65], in the Belle detector at the KEKB [66], and BABAR detector at the PEP-II 
[67, 68]. In addition, at CERN there is a particle accelerator, which is currently 
under construction, namely the Large Hadron Collider (LHC), is expected to 
operate in May 2008 [69]. The collider tunnel in LHC consists of two pipes and 
each pipe contains a proton beam having an energy of 7.0 TeV which travels in 
opposite directions around the ring. This means that in total the collision energy 
will be 14 TeV. Furthermore, the International Linear Collider (ILC), which is 
also under construction, is planned to be completed in the late 2010s [69]. The 
initial collision energy is planned to be 500 GeV and this collision will be between 
electron and positron beams. It is expected that the physics beyond the SM will 
be detected at the LHC and ILC. 

The simplest extension of the SM, the so called the 2HDM, possesses five 
physical Higgs bosons, namely, a charged pair (H^), two neutral CP even scalars 
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[H^ and h^), and a neutral CP odd scalar Observation of these particles in 

the experiments would become a clear signature for the physics beyond the SM. 

Beside this phenomenological hint, the incompatibility of the predicted BRs 
of the LFV interactions with the experimental data obtained, the SM also pos- 
sesses some conceptual problems which motivate us to look physics beyond, such 
as the hierarchy problem between EW and Planck scales. This problem could 
be achieved by introducing the extra dimensions. One possibility is proposed by 
Arkani-Hamed, Savas Dimopoulos and Gia Dvali [20, 21] with n compact extra 
spatial dimensions of large size. In this model, gravity spreads over all the vol- 
ume including the extra dimensions while the matter fields are restricted in 4D 
brane. Another possibility is introduced by Randall and Sundrum [30, 31]. In 
this scenario, the extra dimension is compactified to /Z2 orbifold with two 4D 
brane boundaries. Even if the compactified extra dimensions are so small that 
we perceive the universe as 4D, if there are extra dimensions, fingerprints of them 
sure to exist and such fingerprints are the particles called KK particles which are 
the additional ingredients of the universe with extra dimensions. The masses of 
KK particles are determined by the higher dimensional geometry. In the large 
extra dimensions, for example, the masses of the KK particles are proportional 
to the inverse size of the extra dimension (for two extra dimensions, the size of 
the extra dimension is ~ 1 mm). This tells us that, the current and future ac- 
celerators should be able to discover them. Since in large extra dimensions, the 
only particle that can travel along the extra dimensions is called as graviton and 
thus, it is the only particle that has KK partners. However, the KK partners of 
the graviton interacts as weakly as the graviton itself. Therefore, KK partners 
of graviton would be hard to observe. In warped extra dimensions, however, we 
cannot take just the the inverse size of the extra dimension as the masses of KK 
particles which gives us the Planck scale mass and we know that on the TeV 
brane, nothing much heavier than a TeV can exist. If one calculate the masses 
of KK particles taking into account the warped space-time, the KK particles of 
graviton turn out to have masses of about a TeV and these KK particles, not as 
in the case of large extra dimensions, will interact sixteen orders of magnitude 
larger than the graviton itself in this geometry. This is good news, since depend- 
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ing on the ultimate energy reach in LHC, there is a probabihty of finding the KK 
partners of the graviton. Apart from the graviton, in the Randall Sundrum sce- 
nario, there exists an additional scalar field, that lives in the 5D bulk, such that 
the size of extra dimension is proportional to its vacuum expectation value and 
its fluctuation over this expectation value is called as the radion field. In order 
to avoid the conflict with the equivalence principle, the introduced fleld should 
be massive and to stabilize the distance between the branes, a potential for this 
fleld is proposed by Goldberger and Wise [50] . The radion decays are interesting 
since a considerable information about the scenario under consideration (the RSI 
scenario) can be obtained with the help of accurate measurements. 

In the present work we study the possible LFV decays of the radion fleld r 
and predict the BRs of the LFV r decays r ^ l^l2 in the RSI model. We observe 
that the BRs of the processes we study are at most of the order of 10~^, for the 
small values of radion mass and their sensitivities to decrease with the 
increasing values of m^. On the other hand, the BRs for r ^ (r — > /i^ e^) 

decays are of the order of 10"^^ (10^^^) for the small values of m^. These results 
show that, among these processes, the LFV r ^ /i^ decay would be the most 
appropriate one to measure its BR. 
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Appendix A 



The gauge invariance is an important concept in modern particle theories as it is 
the origin of all of the known four fundamental forces described in Chapter 2. The 
basic method to provide gauge invariance is to ensure that Lagrangian remains 
invariant under certain symmetry transformations which reflect conservation laws 
in nature. By applying these transformations, we end up with conserved physical 
quantities. Since these conserved quantities should not depend on position in 
space-time, theories of particle interactions have to be invariant under local as 
well as global gauge transformations explained below. The transformations could 
be written as 

z/^ ^ t/z/^, (A-1) 
where U is unitarity. So, the symmetry is called U{1) gauge invariance. 



Global Gauge Transformations 

The expression for global gauge transformation (GGT) is 

i/j e^V, (A-2) 

where ^ is a real number. Thus, GGT represents an identical operation at all 
points in space-time and causes a simple shift in the phase of a fermion wave 
function. As a first step we can see how it works in QED. Equation of motion for 
free fermions are obtained from Dirac Lagrangian 

where ip is the fermion spinor. It is clear that such a Lagrangian remains invariant 
under GGT. So, the global phase symmetry is just a statement of the fact that 
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the laws of physics are independent of the choice of phase convention. 



Local Gauge Transformations 



The expression for local gauge transformation (LGT) is 



(A-4) 



where ^ is a function of x = (x, t). Thus, LGT corresponds to choosing a con- 
vention to define the phase of the fermion wavefunction, which is different at 
different points in space-time. Going back to the Dirac Lagrangian defined in eq. 
A-3, one can simply realize that it is not invariant under this more demanding 
symmetry transformation. However, it comes as a pleasant surprise that if we 
introduce another field, A^, a Lagrangian which exhibits local gauge symmetry 
can be obtained. The required field must have infinite range, since there is no 
limit to the distances over which the phase transformations done. Hence, invari- 
ance of Lagrangian requires this new field to be massless. In fact, this field is 
not other than the long-range electromagnetic field: the photon. Therefore, we 
should modify our Lagrangian to make it gauge invariant by replacing the normal 
derivative with the covariant derivative = + iqA^. So, the Lagrangian 
in eq. A-3 reads 

L = ^(27^/^^ - m)ij = ^(27^<9^ - m)i: - qA^i^-f^^^j = Lfree - J^A^. (A-5) 

It so happens that, the photon transforms under LGT as A^i_ ^ A^ + d^9{x) so 
that, the changes in the Lagrangian resulting from the LGT is cancelled out by 
the changes in A^. 

To see the complete picture, we should also add the gauge invariant kinetic 
term to the QED Lagrangian 



LqeD — Lfree — A^ — -F^" 



(A-6) 



89 



where F^j^ is the electromagnetic field tensor defined as 

F^, = d^A, - d,A^. (A-7) 
This Lagrangian expresses the interaction of Dirac fields with the massless photon. 
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Appendix B 



The Einstein Equation in Another Form 

Here we show another form of the Einstein equation which we use in the text. 
By using the equation Gmn = Rain — ^QmnR = i^'^Tmn we obtain Rmn as 

2 1 
Rmn = K, Tmn + -^9mnR, 

and multiplying this equation by g^^^ from left we get 

g^'^'RMN = ^V'Tmn + Ig^'^'gMNR, (B-2) 



1 .„ 2k^ 



which is equal to, 

R = K^T + -dR = ^^T, (B-3) 
2 2-d ^ ' 

where d is the number of dimensions. Substituting this into the eq. B-1 we obtain 

Rmn as 

Rmn — ^ Tmn + -Qmn- — ^-^ 

= k^{Tmn + Qmntt^T) 
2 — a 

= k^Tmn, (B-4) 



with 



where d = 5 in our case. 



Tmn — Tmn — -^OmnT, (B-5) 
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The Derivation of the Linearized Einstein Equa- 
tions 



Now, we consider the metric in eq. 4.110 and first we will derive the explicit 
form of the curvature Rmn and its linearized form. The curvature Rf^i, in terms 
of connection coefficients can be written as 

p _ -pK _ -pK _ -pK -pM I -pK -pM (p, 
^fiu — i- f:,K,u ^ iiu,K ^ nu^ KM ~^ ^ fiM^ uK- \^'^) 

Let US analyze this term by term. The first term is derived as follows 

= ■^{g^^[gLfj.,K + gLK,fi - g^lK,L]},u 

= \{g^''[gp^i,K + gpK,^l - g,,K,p\ + g^^[g5K,t, - gt,K,5]},u 

= \{g^^ [gpf^,^ + gpi,p - ^^s.p] + g^^ [^55,^] } ,^ 



:i + G)-2(l+G)2j,. 
-2e"2^-2^9pFr/^,^) + (1 + Gr^2d,Gil + G)},, 



-4W + fff-^, (B-r) 



and for small fiuctuations we have 



= -Ad^d,F + d^d,G{l-G)-d^GdM^-'^G). (B-8) 



92 



Then, the hnearized form of the first term is obtained as 



The second term is 



515, , = -Ad^d^F + d^d.G. (B-9) 



^t^u,K - {9 Ti^^j^A- 

= -{g^^[gLfi,u + gLu,fj. — gfj.u,L]},K 



7;{g^''[gp^^,y + ^P^,M - ^M^,p] + g^^[-g^^uA},K 



1 
2 

1 1 

1 
2 

+ ^{(l + G)-(e-^-^V).5}.5 

= ^{e^^+^^«1(-2e-^^-^^9.Fr/,,) + (-2e-2^-^^9,Fr/,.) 

- (-2e-2^-2^9pFr/^.)},5 + ^{-2(1 + Gr'd,{A + F)e-2^-2 V},5 

= -rjf^dird^F - Tild^d^^F + r]^Prii,^d^dpF 

_^A-2F , -2{A' + F'f + A' + F" _ 2{A' + F')G'{l^G) ^ 
''^"^ {l + Gf {l + GY ^ 

= -d^d,F-d,dpF + r]p,di:d^F 

^_^^_^^Jh^\ o( A' j_ j_ A// : p// 2(A^ + F')G' 

= -2g,g.F + ry,.DF- ^^^^Jf 

X [-2A'2 - 2F'2 - 4A'F + A" + F" - (B-10) 
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and for small fluctuations we get 



rji,,. = -2d^d,F + 7],,nF 

- e-2V(l-2F)(l-2G') 

X [A" + F" - 2A'^ - 2F'2 - 4A'F' - 2{A' + F')G'{1 - G)].(B-11) 

The linearized form of the second term is obtained as 
dV^^^K = -2d,d.F + v,,OF - e-'^v,. 

[D-12) 

X [F" - AA'F' - 2A'G' - 2FA" + 4FA'^ - 2GA" + AGA% 
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The third term is calculated as 



J- fiu^ KM — 9 ^ Lfj,u9 i SKM 



1 A'L r 1 MS r 1 

= 2^ \.9Lfj,,iy + 9Lu,/i — 9ii,y,L\-9 [9SK,M ^ 9SM,K — 9KM,S\ 

= -^9^^''[9p^L,p + 9pu,^i — 9pu,p\9^^^[9sK,M + 9SM,K — 9km,s\ 

+ -^9^'^^[~9hvA9^^^[9sk,m + 9sm,k — 9km,s] 

= ^9^^[9p^,,u + 9pu,p - 9pu,p]9'''^^[9s(„M + 9sM,i - 9£,m,s\ 

+ -^9^^[—9pu,b]9^^'^[9sb,M + 9sMfi ~ 95M,s\ 

= ^9^''[9pp,u + 9pu,p - 9puA9''^[9Sii,a + 9Sa,l, - 9iia,s\ 
+ \9^^[-9puA9''^[9S^,a + ^5a,5] 

+ ^9^''[9pp,u + - 9pu,p]9^^[9si,b + ^55,d 

+ ^5'^^[-fi'M!-,5]fi'^^[5'55,5 + 9Sb,b - 6-55,5] 
= ^9^''[9pp,u + S'p;.,^ - fl'M^,p]5'"^[5'/3«,a + 9pa,i " S'^Q,/?] 
+ \9''[-9pu,5]9''%^a,5] 

+ ^9^^[9pp,u + - 9pu,p]9^^ [955,i] 

+ ^^7''[-^7m-,5]'7''[^?55,5] 

1 

V 
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= J [-2r]ld,F - 2ri%F + 2r]^,d^F] x [-2ri'^do.F - 8d^F + 2r/f a^F] 
1 16e-^^-'^v,,d,{A + F)d5{A + F) 

+ - [-24d^F - 24d,F + 2r^,^d^F] [2^] 



4^ (1 + G)2 1 + 



Sd^Fd.F - 4r]^,d^Fd^F + 



d^Gd.F d^Fd.G v^^d^Fd^G e-^^~^^r]^,{A' + F')G' 
1 + G 1 + G 1 + G + 



Sd^Fd^F 



1 + G 1 + G (1 + ^)2 



X l-^A' + FT + ^""'^ff ] + V,A-<dFr + '-^], (B-13) 



and for small fluctuations we have 



Cr^V = 8d,Fd,F-d^Gd,Fil-G)-d^Fd,Ga-G) 
- e-^%,{l-2F){l-2G) 
X [-4A'2 - 4F'2 - 8A'F' + {A' + F')G'{1 - G)] 

+ v>.u[-^dF)^ + d^Fd^Gil-G)]. (B-14) 

Then, the linearized form of the third term is obtained as 

K^^u^km) = -e-^\A-8A'F' + A'G' + 8FA'^ + 8GA% (B-15) 
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Finally, we will derive the fourth term in the same manner as 



-pK -pM „KL-p „MS-p 

J- fiM^ uK — 9 ^ LfiAig i SuK 

1 j^'/^ r -1 1 A 'l S r 1 

= 2 ^ [9Lii,M + gLM,^l — giiM,L\ -g [gsv,K + gsK,v — guK,s\ 

= -^g^'^[gp^l,M + gpM^fi — gpM,p]g^^^[gsp,K + gsK,u — gvK,s\ 

+ ~;g [g^M,p ~ g^M, 5]g'''[gsu,K + gsK,u - guK,s] 

= ■^g^^[gpp,M + gpM,p. — gti]\i,p]g'^^^[gsu,^ + gs^,u — gu^,s] 

+ ^g^^[g5M,^^ - gf,Mfi]g^^^[gsy,b + gs^A 



1 



^g^''[gpp,o, + gpa,p - gpa,p]9''^[9sv,i + gs!;,v - g^c 



s 



4 

+ ^g^^"[-gpa,5]9'"'^[gsu,5 + gs5,u] 

+ ^g^^[gpp,5]g^'^[gsu,i + gs^,u - gun, s] 

+ ^g'''[g55,p]9'''^[gsu,5 + gs5A 

= ^g^^[gpp,a + gpa,p - gpajg^'^igpu,^ + gpi,u - gu^A 

+ -/'[-gp.Ag''%puA 

+ \g^''[gppA9''''[-9uiA 

+ h + Gr'il + G)l{l + Gr\l + G)l 



4 
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1 4e-2^-2^r/«r/,595(A + F)95(A + F)^ 



2 



4^ (l + G) 

1 49^^G9^ 

6d,Fd.F-2v,AdFr ' ^V,ue-'^~'^{A' + Fr 



Since we have small fluctuations, we can write 

r^Mr^ = &d,Fd,F-2r^,,{dFf + r^,,,e-^^{l-2F){l-2G) 

X [2A'^ + 2F'^ + AAF']+d^Gd^G{l-2G). (B-17) 

As a result, the linearized form of the fourth term is obtained as 

5{Km^uk) = V,ue-'^[^A'F' - AA"F - AA"G]. (B-18) 

If we add these four terms, we get SR^u, as 

SR^, = -2d^d,F + d^d.G - r]^,aF + e-^^ri^, 

X [F" - 4A'F' - 2A'G' - 2FA" + AFA''^ - 2GA" + AGA'^ 
- 8A'F' + A'G' + 8FA'^ + 8GA'^ + 4A'F' - AA'^F - AA'^G] 
= -2d^d,F + d^dM-r]^,nF + e-'%, 

X [F" -8A'F' - A'G' -2FA" + 8FA''^ -2GA" + 8GA% (B-19) 

Now, we will calculate the curvature i?^5 which can be written in terms of con- 
nection coefficients as 

p _ T^K T^K -pM I pM m on'\ 

-n-AtS — J- imK,5 ii^,K lib'- KM "r ^lM'- 5K- {D-ZV) 
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Let us analyze the -R^s also term by term. The first term is obtained as 



= ■^{g^^[9Lfi,K + gLK,fi — gfiK,L]},5 

= 2^9^^[9p^l,K + QpK,)! — gfiK,p] + 9^^[g5K,p — 9fiK,5\},5 

= ^{9^''[9pf.,i + 9p^,p - 9^^c,p\ + 9^^[955,p]},5 

+ (i + G)-^(i + G)y,5 

_ (_2e-2^-2^9,Fr/^5) + (1 + G)-'2d,G{l + G)},, 

duG , 

5 



{-d,F~Ad,F + d,F+^}, 



5 



Since G = 2F we get 



and taking into account the small fluctuation we obtain 

rf^^5 ^ -49^F' + 29^F'(l-2F)-4F'9^F(l-4F). (B-23) 
Then, the linearized form of the flrst term will be 



Kk,, = -^d,F' + 2d,F' 



2d^F'. (B-24) 
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The second term can be calculated as follows 



^{9'''[9pf.,5]+9''''[955j},K 
1 - ^2A+2Ffpw-2A-2F 



+ ^{(i + G)-^(i + G)y,5 

+ i{(l + G)-22(l+G)9,,G},5 

Substituting G = 2F we have 

W - -^M^ + (rT2F)"(l + 2F)2' (^-26) 
and due to the small fluctuation in F it can be written as 

T%^^ ^ -d^F' + 2d^F'{l-2F)-4F'd^F{l-4F). (B-27) 

We get the linearized form of the second term as 

K^,K = -d,F' + 2d,F' 



d^F'. (B-28) 
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The third term is 



J- ii5^ KM — 9 ^ Lfi59 J- SKM 



X 



1 ji^ij r -| 1 AIS r 1 

= 2 ^ [9Lti,5 + 9L5,ti\ 2 fi' [9sK,M + 9sM,K — 9km,s\ 

= -^9^''[9pnA9^^^[9sK,M + S'SM,^' — 9km,s\ 

+ ^5''^^[5'55,/iJfi'*^'^[fl'5A',M + 9SM,K — 9KM,s] 

= ^9^^[9pfi,5]9^^^[9si,M + 9sM,i - 9iM,s] 

+ ^9'''[955j9 ''[9S5 ,M + fl'SA/,5 — 9bM,s\ 

= ^9^''[9pfi.,5]9''^[9SC,a + 9Sa,^ - 9^a,s] 

+ \9^^[955j9''^[9S5,a + 9Sa,5] 

+ ^9^^[955,p]9^'^[9S5,5 + 9S5,5 - 955,s] 

= ^9^^[9pp,5]9'''^[9l3i,a + 9pa,C " ^^a,/?] 
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^[-2r/«95(A + F)] X [-2r/°a„F - S^^F + 2r/fa^F] 

1 -16a^G95(A + F) 
4^ (1 + G) ^ 

^[-4r/j95(A + F)^] 
1,. d^G d,G ^ 



■[4- 

4^ {1 + G)1 + G' 

Ad^GiA' + F') 



= 4d,F{A' + F') 

d,G{A' + F') G'd.G 
1 + G + 

By using G = 2F we get 



(B-29) 



r^sr^M = ^d,F{A' + F')- +(rT^' ^^"^^^ 

and for small fluctuations we have 

r^r^Af = 49^F(A' + F')-109^F(A' + F')(l-2F) 

+ 4F'9^F(1 -4F). (B-31) 

Then, the third term in the linearized form is obtained as 

'^(rSr^M) = ^A'd.F- 10 A' d,F 

= -QA'd^F. (B-32) 
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The last term is calculated as 



J- fiM^ 5K - 9 ^ LfiMg i S5K 



1 j^Ij r 1 1 A'l S r 1 

= 2^ i9Lti,M + 9LM,ii — 9ij.m,l\-9 [9S5,K + 9sk,5 — 95K,s\ 

= ■^9^''[9pn,M + 9pM,fi — 9pM,p]9^^^[9S5,K + 9SK,5 — 95K,s] 

+ ■^9^^[95M,p — 9pM,^]9^^^[9Sb,K + 9SK,^ — 9bKs\ 

= ^9^''[9pp,M + 9pM,p - 9iiM,p]9^'^^[9sb,i + 9sii,b 

+ -^9^^[9bM,p — 9p,M,5]9^^^[9S5,5 + fi'SS.S — 955,s] 

= ^9^''[9pp,a + 9pa,p - 9pa,p]9''^[9S5,i + ^S^s] 

+ \9^^[-9pa,5]9''^[2gS5,5 - 955,s] 

+ ^^^^[^55,m]^?^^[2^S5,5 -^55,5] 

= ^^^''[^p/.,a + 9pa,p - 9pa,p]9'''^[9Pi,5] 

+ \9''[-9pa,,]9''^[-9,,,p] 

+ \9^'[9pp,5]9'''[955,i] 
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= J[-2r/Ja„F + 2r]id,F - 2rj^ad^F] x [-2r]^d5{A + F)] 

1 -4e-^^-^^r/^^e^-^+^^r/°^a5(A + F)df,G 
4^ (l + G) ^ 

1 -4r/^g5(^ + i^)g€g 4a^Ga5G 
4^ 1 + G (1 + G)2^ 

- 49,F(A+F) ^ — + (rw ^ ^ ^ 

Substituting G = 2F, for small fluctuations, we get 

r^r^ - AaF(A'4^F') <^' + F')d,F AF'd.F 
Wsi. - 4d,F{A+F) 1 + 2F +(1 + 2F)2 

^ 4(9^F(A' + F') - 4(y4' + F')df,F{l - 2F) 

+ AF'd^F{l-AF). (B-34) 

We obtain the linearized form of the fourth term as 

SiT^M^^K) = ^A'd,F-AA'd,F 
= 0. 

If we add these four terms, the linearized form of i?^5 reads 



(B-35) 



SR^5 = -2d^F' -d^F' + QA'd^F + 

= -3df,F' + QA'df,F. (B-36) 

Finally, we will derive the curvature -R55 and, in terms of connection coefficients, 
it reads 

P _ -pK-pM I -pK -pM f-r, Q7\ 

-"-55 — J- 5K,5 '- 55,K ^ 55^ KM "i" 5M^ 5K- ' J 
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Now, we start with the first term in R^^ which reads 



^5^,5 - {9 '^L5k},5 

= -^{g^^[gL5,K + gLK,5 — 95K,l]},5 

= ^{9^^''[9pK,5 - 95kJ + 9^^[955,K + 9bK,b - 5^5^,5] },5 

= ^{/1^?pC,5] +/'b5,5]},5 

= ^{e^^^^V''(e-^^-^^p,),5 + (1 + G)-(l + G)l], 

- + n + (B-38) 

Since G = 2F we get 

rf,, = -4(^» + n + ^-(^^. (B-39) 
For small fluctuations we have 

^ -4(A" + F") + 2F"(1 - 2F) - 4F'2(1 - 4F). (B-40) 
Then, one can simply write the linearized form of the first term as 

5r^5 = -AF" + 2F" 



-2F". (B-41) 
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The second term is 



^{9^'[-955,p] + ^^^[2^55,5 - ^55,5]},K 
^{/l-^755,p]},5 + ^{/'(^?55,5)},5 

e^^+^Fl2d^Fd^Gil + G)+ d^d^Gil + G) + d^Gd^G] 
G" G'2 



1 + G (1 + G)2- 
Substituting G = 2F we get 



(B-42) 



Tg^^ = e2^+'^[4(9F)2(l + 2F) + 2nF(l + 2F)+4(aF)2] 

+ TT^- {i+2Fr ^^'^^^ 

and for small fluctuations we have 

Tgj^ ^ e2^(l + 2F)[4(9F)2(1 + 2F) +2nF(l + 2F) +4(9F)2] 

+ 2F"(1 - 2F) - 4F'2(1 - 4F). (B-44) 

Then, the linearized form of the second term is obtained as 



Srg K = 2e^^OF + 2F". (B-45) 
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The third term in the curvature -R55 is 



J- 55J- KM — 9 ^ L559 J- SKM 



1 r 1 MS r 1 

[9L5,5 + 9L5,5 — 955,l\ -9 [9SK,M + 9SM,K — 9KM,s\ 



\9'"'[-9^^j9'''[9sK,. 
1 



L J^'^^yjo ^j^s.,M + 9SM,K — 9KM,S\ 



+ ■^9^^['^955,5 — 955,5]9'^^[9SK,M + 9SM,K — 9KM,s] 
^9^^[-955,p]9^^^[9Si,M + 9SM,^ - 9iM,s] 
-^9^^[955,5]9^^^[9S5,M + 9SM,5 — 95M,s] 



^9^^[-955,p]9"^[9Si,a + 9Sa,i - 9ia,s] 
1 



+ ^9^^ [5'55,5]fl'"^ [9S5,a + 9Sa,5] 

+ \9^''[-955,p]9^^[9Si,5+9S5,i] 

+ ^9^^ [5'55,5]fl'^"^ [9S5,5 + 9S5,5 " 955,s] 

= \9^^[-955,p]9"'^[9l3i,a + 9Pa,i " 9^cA 

+ -/'[9^^,^K%PaA 

+ ^(1 + G)-2(1 + G')25(1 + G)-2(1 + G)25 
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1 jg2A+2F2(l + G)d^G] X [-27]P^F - 8d^F + 2r]^dpF] 



1 -16G'd,{A + F) 
4^ (l + G) ^ 

4^ ^ ^ l + G 

1 . . G' d^G' ^ 



4^ 1 + + 

^ l + G) 



+ e^^+^^a^Ga,G + ^^^, (B-46) 
and by using G = 2F, the third term is obtained to be equal to 

+ 4e^^^-^{dFr + j^^^^, (B-47) 

and, considering the small fluctuation, we get 

rfsr^M = -86^^(1 + 2F)(1 + 2F)(9F)2-8(A' + F')F'(1-2F) 

+ 46^^(1 + 2F)(9F)2 + 4F'2(1-4F). (B-48) 

Then, the linearized form of the third term is obtained as 

Sirg^KM) = -8A'F'. (B-49) 
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Finally, the fourth term is derived as follows: 



J- 5mJ- 5K — 9 ^ L5M9 i S5K 



1 A'L r 1 A-IS r 1 

-^9 [9l5,m ^ 9lm,5 — 95M,l\-9 [9s5,k + 9sk,5 — 95K,s\ 

■^9^''[9pM,5 — 95M,p]9^^^[9S5,K + 9SK,5 — 95K,s] 
^9^^[955m]9^^^[9S5,K + 9SK,5 - 95KS] 

1 

4 

^9^^[955,m]9'^^^[9S5,5 + 9S5,5 - 955,s] 



-;9^^[9pM,b - 9bM,p]9^'^[9s^,^ + 9siA 



+ \9^'[9^5A9^^[9su + 9siA 
= -/'[9p.A9'''[9kA 



1 



+ ^(1 + G)-^(1 + G)25(1 + G')-^(1+G)^5 
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+ 
+ 



l[-27^id,{A + F)][-2r^^d,{A + F)] 

1 4e^^+2V^(l + G)daGdf3G 
4^ 1 + G ^ 

1 Ae'^^+^^d^Gd^Gjl + 
4^ 1 + G 

1.. G'2 , 



+ 7[4 



4^ (1 + G)2J 

= 4(A' + F')' + 2e2^+2^(9G)2 + — (B-50) 

(1 + Gj 

Substituting G = 2F we obtain 

rfMrf^ = 4(A' + F')^ + 8e2^+2^(9F)^ + ^^i^, (B-51) 
and, with small fluctuations, we have 

rf^r^ ^ 4(v4' + F')' + 8e2^(l + 2F)(9F)2 + 4F'2(l-4F). (B-52) 
Then, the linearized form of the fourth term is obtained as 

SirfM^fK) = 8A'F'. (B-53) 

As a final step we will add these four terms so that we get the linearized form of 

-R55 as 

= -4F" - 2e2^nF + IQA'F'. (B-54) 

At this stage, we will derive the linearized form of the source term. We start with 
the part of the action (see eq. 4.78) including the source term and taking ^/g55 
as ^/g55 = 1 + 2F. The metric variation on this action gives 



(B-55) 
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Since T = QmnT we get 

i 

Using the form of Einstein equation given in eq. B-1, and the eqs. B-55 and B-56 
to get 

T^y = T^u - -gf,uT, (B-57) 

we obtain 

(B-58) 

where 

0(x,y) = (/)o(y) + </?(x,y), 

\/(</.) = \/(</.o) + <^y'(0o), (B-59) 
Ai(0) = Ai(</)o) + V5Ai(0o), 

for small fluctuations. Substituting the equations in eq. B-59 into Ty,y we have 

1 ^ (B-60) 
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The linearized form of the source term is obtained as 

5%, =ie-2^r/^,[<^V'(0o) - 2V{<t)^)F] 



Now, let us derive the equation for the source term T„5 by using 



~ _ 1 

^At5 — T'^S — -^9fi5T. 



Since (7^5 = we get 



Tfj,5 — T^5. 



Using the eq. B-55 we obtain T^s as 



T^5 = -^d^c^d^^ 



'^[dMoiy) + v^ix, y))d5{My) + vi^^ y))]- 



Then, the linearized form of T^s reads 



Finally, T55 can be obtained by using the eq. B-55 
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Then, the source term T55 becomes 

^55 = ^55 — gfl'SS^ 



1 1 5 (B-67) 

+ -(l + 2F)2[-(a0)2-|l^(0)] 

-hl + 2F)Y,Hmy-y^)■ 



Making the simphfications and substituting the eqs. B-59 we get 

T55 = -^(1 + 4F + 4F2) [V{<P^) + vV\<P,)] 

- ]^db[<Po{y) + v{x,y)%[(t)o{y) + ^{x,y)] (B-68) 

- |(1 + 2F) 5^[A,(0o) + vK{<P,my - yi). 



The hnearized form is obtained as 



5f,, = -U^{<P,)F-^vV\<P,) 
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Appendix C 



In this Appendix we present the formulation for the spin connection which we 
use in our calculations (see [70]). A natural basis for the tangent space Tp at 
a point p is given by the partial derivatives with respect to the coordinates at 
that point, e(^) = d^. Similarly, a basis for the cotangent space T* is given by 
of the coordinate functions, = dx'^. Let us imagine that at each point in the 
manifold there exists a set of orthonormal basis vectors e(a)^. If the canonical 
form of the metric is written r]ab, the inner product of our basis vectors should be 

e{a) ■ 6(6) =rjab. (C-1) 

Thus, in Lorentzian space rjat represents the Minkowski metric, while in a space 
with positive definite metric it represents Euclidean metric. We can express our 
old basis vectors of tangent space in terms of the new ones as 

where the components are called as vielbeins form an invertible n x n matrix. 
Their inverse is denoted by switching the indices to obtain e^, which satisfy 

e^^et = 6i: ■ ey,=5l (C-3) 

Multiplying the eq. C-2 with from left 

e(a)=e^e(^), (C-4) 

^indexed in Latin letter rather than Greek, to remind us that they are not related to any 
coordinate system. 



114 



Then, in terms of the inverse vielbeins the eq. C-1 becomes 



Vab = e(a) ■ 6(6) 

= <e(^) ■ e^e(,) (C-5) 



We see that the components of the metric tensor in the orthonormal basis are 
just those of flat metric, rjab- Multiplying this with e^e^ we get 



9,.u = e^^e%b- (C-6) 

Thus, vielbeins are squareroot of the metric. Similarly in T*, 9^"^^ are the or- 
thonormal basis of one-forms. Choosing 

^('^)e(,) = S^, = elel (C-7) 

it is an immediate consequence that the orthonormal one-forms are related to the 
cooridinate-based ^(^) by 

Similarly, 

ka) = e%). (C-9) 

Any vector V written in the coordinate basis as ^ = V^e^ can be expressed as 
in terms of its orthonormal basis as \^ = V^Ca- Since V^e^ = V^Ca we obtain a 
relation between the sets of components as 

So the vielbeins allow us to pass from Latin to Greek indices and back. In the 
same manner multi index tensor Vj," can be written as 
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Nice property of tensors is that, we can go on to refer multi index tensors in terms 
of mixed components. 

Now, we have a set of basis vectors Ca and Oca which are non-coordinate 
bases. Thus, they can be changed independently of the coordinates provided the 
orthonormahty property defined in eq. C-1 is preserved. In Euchdean signature 
metric the transformations that preserve orthonormahty condition are orthogonal 
transformations whereas in Lorentz signature metric they are Lorentz transfor- 
mations. We therefore consider changes of basis of the form 

Ca ^ ea> = A/(x)e„ (C-12) 

where A", (x) represent position dependent transformations at each point in space 
which leave the canonical form of metric unaltered such that 

K,,\x)\,\x)7]ab = Va'b'. (C-13) 

In flat space, we call these matrices inverse Lorentz transformations. We also 
have ordinary Lorentz transformations, A^'(x) to transform one-forms. So, we 
now have freedom to perform a Lorentz transformation at every point in space. 
These are called local Lorentz transformations (LLT). 

The covariant derivative of a tensor is given by its partial derivative plus 
the connection terms, one for each index. The connection terms in our ordinary 
formalism involve the tensor and connection coefficients F^^^ whereas in non- 
coordinate basis connection coefficients are replaced by spin connection, denoted 
by Each Latin index gets a factor of the spin connection as 

v,x\ = d,x\ + w; , - w; ,x\. (c-i4) 

Now let us try to find a relation between the spin connection, the vielbeins and 
connection coefficients using the property that a tensor should be independent of 
the way it is written. For simplicity, we will consider the covariant derivative of 
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a vector X. Its covariant derivative in a purely coordinate basis is given by 



VX = (V^X'^)c/x^ (g) d, 
and in a mixed basis it is written as 



(C-15) 



VX = {V^X')dx^ ® e(a) 

= {d^X'^ + w;,X'^)dx^®e^a) 

= {d,{elX'^)+w;AX^)dx^®ha) 

= {eld.X-' + X-'d.e: + w; Ax')dx^ ® e^d^ (C-16) 
= elield^X" + X'd^et + ,e\x'')dx^ ® 8, 
= {8ld^X^ + elX^d^el + e>/ A^^^^x^ ® 
= {d^X- + elX'^d^el + e>/ ,e^X^)rfx'^ ® d^. 

Let cr ^ z/ and u ^ X 

VX = (a^X'^ + e^(9^e^)X^ + e^e^w;/ ,X")t/x^ ® d,. (C-17) 
Comparing the eqs. C-15 and C-17 one can conclude that 

T'',, = e':{d,el)X^ + e':e\w;,. (C-18) 
Multiplying this with e^e" we get 

w;, = elelTl^-eld,el. (C-19) 
Now let us look the covariant derivative of a vielbein 

V,el = d,el + w; ,el - Y^e^. (C-20) 
Substituting the equation for V^^y^ into the covariant derivative of the vielbein we 
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(C-21) 
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